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ABSTRACT 

Although the spectrum of a prototypical early-type galaxy is assumed to lack emission lines, a 
substantial fraction (likely as high as 30%) of nearby red sequence galaxy spectra contain emission 
lines with line ratios characteristic of low ionization nuclear emission-line regions (LINERs). We 
use spectra of ~ 6000 galaxies from the Sloan Digital Sky Survey (SDSS) in a narrow redshift slice 
(0.06 < z < 0.08) to compare the stellar populations of red sequence galaxies with and without 
LINER-like emission. The spectra are binned by internal velocity dispersion (a) and by emission 
properties to produce high S/N stacked spectra. The recent stellar population models of R. Schiavon 
(2006) make it possible to measure ages, [Fe/H], and individual elemental abundance ratios [Mg/Fe], 
[C/Fe], [N/Fe], and [Ca/Fe] for each of the stacked spectra. We find that red sequence galaxies 
with strong LINER-like emission are systematically 2-3.5 Gyr (10-40%) younger than their emission- 
free counterparts at the same a. This suggests a connection between the mechanism powering the 
emission (whether AGN, post-AGB stars, shocks, or cooling flows) and more recent star formation 
in the galaxy. We find that mean stellar age and [Fe/H] increase with a for all galaxies. Elemental 
abundance [Mg/Fe] increases modestly with a in agreement with previous results, and [C/Fe] and 
[N/Fe] increase more strongly with a than does [Mg/Fe]. [Ca/Fe] appears to be roughly solar for all 
galaxies. At fixed a galaxies with fainter r-band luminosities have lower [Fe/H] and older ages but 
similar abundance ratios compared to brighter galaxies. 
Subject headings: galaxies: abundances, galaxies: elliptical and lenticular 



1. INTRODUCTION 

An important recent result in the field of galaxy evo- 
lution is the build-up of the red sequence between 2 = 1 
and the present. The red sequence luminosity func- 
tion determined by the COMBO-17 survey (Bell et al. 
2004) shows an accumulation of massive red galaxies 
since z = 1. This result was confirmed by Faber et al. 
(2006) using combined data from COMBO-17 and the 
DEEP2 redshift survey. Faber et al. (2006) contrasted 
the evolving luminosity function of red galaxies with the 
luminosity function of massive blue galaxies, which re- 
mains roughly constant over the same epoch. With so 
much transition in the relatively recent past, there should 
be visible evidence of recent red-galaxy creation in the 
local universe; that is, nearby galaxy samples might con- 
tain galaxies that are currently evolving onto the red 
sequence, as well as recent red sequence arrivals. 

Stellar population models predict that once a galaxy 
has ceased to form stars it will move to the red sequence 
rapidly. A galaxy whose star formation is quenched sud- 
denly after a history of continuous constant star forma- 
tion will arrive on the red sequence within ~ 400 Myr, 
while a galaxy that experiences a substantial burst of 
star formation before quenching will take only ~ 1 Gyr 
to become red (Harker et al. 2006). The red sequence 
may therefore include galaxies which contain a fractional 
population of stars with young ages. Under the possi- 
bly more reasonable assumption of star formation that 
declines with time, the fraction of young stars will be 
smaller, will contribute less to the total galaxy light, and 
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presumably the galaxy will move to the red sequence 
more rapidly. 

Recently, Yi et al. (2005), Kaviraj et al. (2006), 
and Schawinski et al. (2006) have used near-ultraviolet 
(NUV) photometry from the GALEX Medium Imag- 
ing Survey in conjunction with Sloan Digital Sky Sur- 
vey (SDSS) optical photometry to examine the colors of 
early-type galaxies. They show that morphologically se- 
lected early-type galaxies have uniformly red optical col- 
ors but exhibit a wide range of NUV-optical colors. Ac- 
cording to their analysis, the observed blue NUV-optical 
colors in some early-type galaxies are indicative of low 
levels of residual star formation; they use Bruzual & 
Chariot (2003) stellar population models to demonstrate 
that blue ./VC/U-optical colors are sensitive to very low lev- 
els of recent star formation and determine that at least 
30% of red sequence galaxies have 1-3% of their mass 
in young (age < 1 Gyr) stars, with typical ages for the 
young component of 300-500 Myr. 

Although the NUV-optical colors of these morpholog- 
ically selected galaxies show signs of recent star forma- 
tion, the galaxies in their sample lack optical emission 
lines, almost without exception. Kaviraj et al. (2006) 
show that their morphological selection criterion removes 
virtually all galaxies with emission lines characteristic 
of Hn regions. In addition, they exclude all galaxies 
with emission lines characteristic of active galactic nu- 
clei (AGN) and low ionization nuclear emission regions 
(LINERs) , using the criterion established by Kauffmann 
et al. (2003). This is to remove from the sample galaxies 
that might have substantial ultraviolet (UV) light gen- 
erated by a source other than the stars. However, by 
excluding all galaxies with AGN-like emission lines, they 
may have removed many red sequence galaxies from their 
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sample which are not forming stars. These galaxies could 
include many of the most recent arrivals on the red se- 
quence. 

Kauffmann et al. (2003) show that the host galaxies of 
luminous ([Oiii]A5007 emission line luminosities in excess 
of 10 7 L Q ) AGN appear to be young or post-starburst sys- 
tems, while those of low-luminosity active nuclei (a cat- 
egory dominated by LINERs in their sample) have mor- 
phologies and stellar populations that are comparable to 
those of early-type galaxies with no emission. Recently, 
Yan et al. (2006) have shown that up to 30% of red se- 
quence galaxies in the SDSS have emission lines, with 
line ratios characteristic of LINERs. Another ~ 50% of 
the red sequence consists of "quiescent" galaxies with no 
detectable emission, while the remaining ~ 20% are di- 
vided among Seyferts, dusty star-forming galaxies, and 
Transition Objects, which likely contain both Seyfert 
and Hll region emission (Yan et al. 2006). A signifi- 
cant fraction of the red sequence population therefore 
consists of galaxies with optical line emission, although 
these galaxies are not necessarily forming stars. These 
galaxies are often systematically excluded from stellar 
population analysis of red sequence or early-type galax- 
ies, even though they are a large class of objects and 
may harbor clues to the evolution of galaxies onto the 
red sequence. 

This paper presents a stellar population analysis of 
red sequence galaxies, including both quiescent galax- 
ies with no emission and those with LINER-like emis- 
sion. The fact that galaxies with LINER-like emission 
lines make up 30% of the red sequence, in addition to 
having morphologies and stellar populations that appear 
to be consistent with the larger population of early-type 
galaxies generally, justifies a stellar population analysis 
performed in parallel with quiescent red sequence galax- 
ies. Also, if LINERs are genuinely active nuclei and if 
their host galaxies form the low luminosity end of the 
AGN host galaxy sequence described in Kauffmann et al. 
(2003), these galaxies with LINER-like emission may be 
candidates for "post post-starburst" systems that have 
recently stopped forming stars and moved onto the red 
sequence. Signs of such recent quenching could then be 
detected through a careful differential stellar population 
analysis of LINERs versus non-LINERs. 

There is still no consensus on the physical mechanism 
that powers LINER emission. Possibilities include pho- 
toionization by a central AGN accretion disk as in higher 
ionization Seyferts (Ferland & Netzer 1983, Halpcrn & 
Steiner 1983), excitation through shocks (as in the mod- 
els of Dopita & Sutherland 1995), accretion of warm gas 
in cooling flows (e.g. Heckman 1981), or photoionization 
by post asymptotic giant branch (AGB) stars in either 
old (Binctte et al. 1994) or young (Taniguchi, Shioya, & 
Murayama 2000) stellar populations. This topic is cov- 
ered by the Filippcnko (2003) review, which includes a 
nice summary of these mechanisms. 

There are clearly LINERs whose emission is strongly 
dominated by low luminosity cental AGN, identified by 
either a central non-stellar source (an X-ray or UV point 
source, or a compact radio core) or by a signature of AGN 
accretion (broad line emission or even double-peaked 
broad Ha emission) . Examples of each of these are refer- 
enced in Filippenko (2003). However, there is also sub- 
stantial evidence for ionized gas in early type galaxies 



that extends multiple kiloparsecs from the galaxy center, 
as shown in narrow-band imaging of early-type galax- 
ies (e.g., Goudfrooij et al. 1994, Macchetto et al. 1996) 
and more recently in the SAURON integral field unit 
spectroscopy of elliptical and lenticular galaxies (Sarzi 
et al. 2006). These studies imply that a large fraction 
of luminous elliptical galaxies contain distributed ionized 
gas with optical emission line ratios that are character- 
istic of LINERs. We must be careful with terminology, 
since the original definition of the term LINER (Heckman 
1980) required that the emission be nuclear ("low ioniza- 
tion nuclear emission-line region"). These galaxies are 
therefore not nuclear LINERs, but they do exhibit emis- 
sion with LINER-like ratios. We shall therefore refer to 
them as "LINER-like" galaxies, without any assumption 
about the spatial distribution of the ionized gas within 
the galaxy or the ionization mechanism. From SDSS data 
alone, whose spectral fibers cover a substantial portion of 
the galaxy, it is not possible to decide whether any given 
galaxy is a classic AGN LINER, an extended LINER-like 
galaxy, or some combination of the two. 

The fact that LINER-like emission is common in red 
sequence galaxies prompts many questions: is it an evo- 
lutionary phase that all early-type galaxies pass through? 
Do galaxies pass through this stage more than once? Per- 
haps many times? Is it a cyclical process? Is the gas ex- 
ternal in origin, as might be suggested by the kinematic 
decoupling of gas and stars seen by Caon, Macchetto, 
& Pastoriza (2000) and Sarzi et al. (2006), or is it pro- 
duced by mass-loss from stars within the galaxy? What 
is the relation, if any, between the central AGN LINER 
emission and the more extended LINER-like emission? 
In this paper, we seek to address some of these questions 
by examining the stellar populations of the host galaxies 
to determine whether or not the existence of LINER-like 
emission in a red sequence galaxy is related to its star 
formation history. 

Past stellar population analysis of early-type or red se- 
quence galaxies has tended to either exclude all galaxies 
with emission (e.g., Nelan et al. 2005 and Bernardi et 
al. 2003b) or to include galaxies with weak LINER-like 
emission in the main sample. In the latter case, authors 
have made a correction for Balmer line emission infill 
using an empirical correction based on the strength of 
the [Om]A5007 (e.g., Gonzalez 1993, Trager et al. 2000a, 
Rampazzo et al. 2005) or [On]A3727 (e.g., Schiavon et 
al. 2006) emission lines, or else restrict themselves to us- 
ing higher-order Balmer lines such as H7 and H<5, where 
contamination from emission is less. 

In this work, we analyze galaxies with and without 
LINER-like emission in parallel, in order to contrast their 
stellar population properties. LINER-like nebulae pro- 
duce only modest quantities of Balmer emission, and that 
emission is well-characterized (due to the selection pro- 
cess described in §2.2). This allows us to make accurate 
Balmer infill corrections for multiple Balmer lines and 
to compare results. In order to obtain very high signal 
to noise (S/N) spectra for the stellar population analy- 
sis, we stack several hundred similar galaxies to produce 
highly accurate determinations of the mean stellar pop- 
ulation properties of the galaxies that go into each com- 
posite. We use the new single stellar population (SSP) 
models of Schiavon (2007), which allow the determina- 
tion of not only mean ages and total metallicities, but 
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of several individual elemental abundances as well; iron, 
magnesium, carbon, nitrogen, and calcium are all sepa- 
rately determined. Section 2 presents the data sample, 
the criteria for grouping similar galaxies used to con- 
struct the composite spectra, and the procedure used 
to correct for Balmer emission infill in the LINER-likc 
galaxies. Section 3 gives results of our line strength mea- 
surements as a function of galaxies properties. Section 
4 describes the stellar population modeling process and 
important caveats, then §5 uses SSP models to interpret 
the line strength measurements in terms of stellar pop- 
ulation properties. In §6, we discuss the implications 
of our results on the possible source of the LINER-likc 
emission in red sequence galaxies, and §7 summarizes our 
conclusions. 

2. DATA 

This paper presents an analysis of red sequence galaxy 
spectra from the Sloan Digital Sky Survey (SDSS), a 
massive imaging and spectroscopic survey conducted at 
the Apache Point Observatory with a dedicated 2.5 m 
telescope. The SDSS has imaged more than one quar- 
ter of the sky in five filter bands: u, g, r, i, and z 
(Fukugita et al. 1996, Stoughton et al. 2002). The survey 
has also obtained spectroscopic observations of over one 
million galaxies and quasars using a purpose-built dual 
fiber spectrograph, with resolution R = A/AA s=s 1800, 
wavelength coverage A = 3800-9200A, and fiber diame- 
ter d — 3". The data presented in this paper are taken 
as part of the SDSS main galaxy sample, which includes 
over 500,000 galaxies with redshifts 0.00 < z < 0.30. 
The main galaxy sample is selected to include galaxies 
brighter than r < 17.77. 

The following imaging and spectroscopic parameters 
for each galaxy were taken from the NYU Value Added 
Catalog (Blanton et al. 2004, Adclman-McCarthy et al. 
2005), Data Release 4 (DR4): rcdshift z, velocity dis- 
persion cfib as measured in the SDSS fiber spectrum, 
ugriz fluxes measured in the SDSS flux units "mag- 
gies" (see Blanton et al. 2004), and median S/N per A 
for each spectrum. The ugriz fluxes were converted to 
absolute magnitudes assuming standard ACDM cosmol- 
ogy (CIa = 0.7, TO = 0.3), with a Hubble constant of 
Hq = 70 km s _1 Mpc -1 . The magnitudes were then 
K-corrected to z = 0.1 using Mike Blanton's kcorrect 
code, v. 3_2 (Blanton ct al. 2003a). Wc chose to In- 
correct to z = 0.1 both to enable easy comparison with 
colors and magnitudes for SDSS galaxies in Blanton et 
al. (2003b) and because we limited our galaxy sample to 
galaxies with 0.06 < z < 0.08 (see §2.1); the K-correction 
to z = 0.1 is small and therefore any errors due to un- 
certainty in the K-correction are minimized. Through- 
out this paper, all magnitudes and colors are quoted at 
z = 0.1, unless otherwise noted. 

These values were supplemented with the following 
measurements from the SDSS DR4 Catalog Archive 
Server: the best-fit deVaucouleurs radius (r e ) and axis 
ratio {b/a). Our analysis also uses emission line equiva- 
lent widths (EWs) for Ha, H/3, [On]A3727, [Om]A5007, 
[Oi]A6300, and [Nn]A6583 measured in SDSS DR4 galax- 
ies by Yan et al. (2006), with the zeropoint corrections 
defined in their Tabic B3. Before making their emission 
line measurements, Yan et al. fit and subtract off a stellar 
absorption line template spectrum to obtain good mea- 



surements of even those emission lines that fall on top of 
stellar absorption features. The interested reader should 
consult their paper for further details on this process. 

2.1. Sample Selection 

The sample is chosen to match the selection criteria of 
the Yan et al. (2006) magnitude limited sample; objects 
are from the SDSS main galaxy sample and are spectro- 
scopically confirmed to be galaxies. We then limit the 
sample to a narrow range in redshift, 0.06 < z < 0.08, 
making this magnitude-limited sample roughly volume- 
limited and minimizing the impact of evolution on this 
analysis. The 1.5" SDSS spectral fiber radius corre- 
sponds to 1.7-2.3 kpc in this redshift range. 

The aim of this work is to study red sequence galax- 
ies with no active star formation. We use two selection 
criteria to select galaxies which are not currently form- 
ing stars: a color cut (described below) and the emis- 
sion line properties of the galaxies. Yan ct al. (2006) 
demonstrated that the distribution of [On] to Ha emis- 
sion line EW ratios is bimodal such that some emission 
line galaxies have strong Ha emission but relatively weak 
[On] emission while others have strong [On] and weak 
Ha. This is best seen in a plot of log EW(Ha) vs. log 
EWQOii]), as in Figure la (see also Yan et al. 2006, Fig- 
ure 2). 

They showed that the left peak with high EW([On]) 
andlowEW(Ha) (they call these "High- [On] /Ha" galax- 
ies) is dominated by red sequence galaxies, whereas the 
right peak with larger EW(Ha) ("Low- [On] /Ha" galax- 
ies) is dominated by blue galaxies. This led them to 
explore the possibility that the High- [On] /Ha galaxies 
are star-forming systems while the Low- [Oil] /Ha galax- 
ies contain emission that is not powered by star forma- 
tion. They confirm this hypothesis through the use of 
line ratio classification diagrams first proposed by Bald- 
win, Phillips, & Terlevich (1981) and commonly known 
as "BPT diagrams" . 

Following Yan et al. (2006), we define red sequence 
galaxies as those with 

01 (.9 - r) > -0.025( 0A M r - 5 log h) + 0.42, (1) 

here with h = 0.70. This includes 22,501 galaxies in the 
chosen redshift range. This color cut is shown in Figure 
lc as the dashed line. The color cut has been chosen 
to be conservative and excludes galaxies that lie in the 
"valley" between the red sequence and the bulk of blue 
galaxies. On the red sequence, we use Ha and [On] EWs 
to characterize galaxies according to their emission line 
properties. We use the Yan et al. (2006) relation 

EWQOii]) > 5EW(Ha) - 7, (2) 

to separate High- [On] /Ha from Low- [Oil] /Ha galaxies. 
This dividing line is shown as the solid line in Figures la 
and lb. Figure la shows the same information as Figure 
lb, but with linear axes. The latter format does not 
highlight the bimodality as clearly but is a more natural 
way to examine the cut used to define the High- [On] /Ha 
sample. 

To confirm that the High- [On] /Ha objects are not star 
forming galaxies, Figure 2 shows BPT diagrams of emis- 
sion line galaxies in the redshift range 0.06 < z < 0.08. 
Red points show the red sequence High- [On] /Ha objects 
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as defined above; black points show all other emission- 
line galaxies. Only galaxies with a > 3c detection in 
all lines are shown on each plot. Panel (b) contains 
fewer objects than (a) because [Oi]A6300 is a weak line 
and undetected in 2/3 of the High- [On] /Ha and 1/2 
of all other galaxies shown in (a). The solid lines in 
each panel show the standard classification scheme di- 
viding Scyfcrts, LINERs, and star forming galaxies. The 
dotted lines show the Kewley et al. (2001) division be- 
tween star forming galaxies and active nuclei; the dashed 
line in (a) shows the division used by Kauffmann et al. 
(2003). The dash-dot line in (b) shows the Ho, Filip- 
penko, & Sargent (1997) division between LINERs and 
transition objects (TOs). Based on the Kauffmann et 
al. (2003) division between star forming galaxies and ac- 
tive nuclei, the standard division between LINERs and 
Seyferts, and the Ho, Filippcnko, & Sargent (1997) divi- 
sion between LINERs and TOs, the red sequence High- 
[On]/Ha galaxies consist of 88.4% LINERs, 8.6% TOs, 
1.9% Seyferts, and 1.1% star forming galaxies. Most of 
the High- [On] /Ha objects classified as TOs lie very near 
the division between LINERs and TOs, well within the 
typical errors, and are not convincingly a separate class 
of objects. There arc High- [On] /Ha galaxies that are 
not shown here because their H/3, [Om]A5007, [Oi]A6300, 
and/or [Nn]A6583 emission lines are weaker than 3cr de- 
tections, but they are assumed to possess the same un- 
derlying emission line properties as those shown. 

The analysis in this paper uses two samples of galaxies 
classified by their [On] and Ha emission lines: galax- 
ies with neither Ha nor [Oil] detected (at the 2<r level) 
and High- [On] /Ha galaxies with both lines detected. We 
shall hereafter refer to these as "quiescent" and "LINER- 
like" galaxies respectively. The color distributions of 
these samples are similar and concentrated on the red 
sequence. Figure Id shows CMDs of the quiescent and 
LINER-like galaxies, with the color cut defining the red 
sequence overplotted as the dashed line (as in Figure lc). 
While the colors of the LINER-like and quiescent galaxies 
are similar, the luminosity distribution is clearly differ- 
ent. This is partially a selection effect; the LINER-like 
galaxies are required to have both [On] and Ha detected 
at the 2a level, thus fainter galaxies with lower S/N spec- 
tra are less likely to meet the selection criteria. Taking 
this effect into account, there still appears to be a de- 
cline in the fraction of red sequence LINER-like galaxies 
at faint magnitudes. This will be explored in more detail 
in future work. Only galaxies redder than the dashed 
line of Figures lc-d are included in our sample. Spectra 
of the quiescent and LINER-like red sequence galaxies 
were downloaded from the SDSS Data Archive Server 
(http: / / das.sdss.org) . 

As defined, our sample of galaxies contains overlap 
with the sample of AGNs and their hosts described by 
Kauffmann et al. (2003). They select galaxies on the ba- 
sis of their log([Nn]A6583/Ha) and log([Om]A3727/H/3) 
ratios, including galaxies that are classified as Seyferts 
by the criteria of Ho, Filippcnko, & Sargent (1997), as 
well as LINERs and some Transition Objects. They 
claim that all of these objects harbor AGN of some 
kind, and divide their sample into low-luminosity and 
high-luminosity AGN, with the dividing line at log 
£[OlH]/£© = 7.0. Careful study of their Figure 3 re- 



veals that at least 85% of the objects that they classify 
as LINERs fall into their low-luminosity AGN sample, 
while only ~ 20% of their Seyferts fall into this sample. 
Comparisons to the results of Kauffmann et al. (2003) 
are therefore most appropriate between our sample and 
their "low-luminosity AGN" sample. However, they ap- 
ply no color criterion to their selection of AGN hosts. Of 
the galaxies in our Figure 2 which meet the Kauffmann 
et al. (2003) definition of AGN, 55% are bluer than our 
red sequence cut-off and 37% meet the more conservative 
color criterion for blue galaxies in Yan et al. (2006). The 
Kauffmann et al. (2003) sample clearly includes many 
more young (and possibly actively star forming) galaxies 
than our sample. Their analysis of other properties of the 
host galaxies suggest that the hosts of strong AGN are 
young galaxies but that the hosts of their low-luminosity 
AGN sample are similar to early-type galaxies, so many 
or most of these blue galaxies may be in their high- 
luminosity AGN sample. However, their low-luminosity 
AGN sample is still likely to contain a significant fraction 
of galaxies bluer than any present in our analysis. 

2.2. Emission Infill Correction 

Our analysis uses the strength of the stellar H/3, H7, 
and H<5 absorption features, in conjunction with metal 
lines, to measure SSP ages for the galaxies in our sam- 
ple. In the LINER-like galaxies, stellar Balmer line ab- 
sorption will be partially or completely filled in by emis- 
sion from the ionized gas within the galaxy. This infill 
severely impacts stellar population analysis, and a cor- 
rection must be made. If we use the Ha and H/3 emission 
line strengths measured by Yan et al. (2006) to make an 
emission infill correction, we merely recover the Balmer 
absorption line strengths of their subtracted absorption 
line spectral template, and further stellar population 
modeling will be meaningless. Gonzalez (1993) defined 
an H/3 infill correction based on the EW of [Om]A5007 
emission, which has been used (sometimes in modified 
form) by many other groups. Trager et al. (2000a) use 
AH/3 = 0.7 EWQOin]) which, assuming the continuum 
level is similar at H/3 and [Oni], is a good match to 
the observed [Oiii]/H/3 ratios for LINERs shown in Fig- 
ure 2. However, [Oin]A5007 is weaker in LINERs than 
[Oli]A3727 due to their lower ionization. For the SDSS 
spectra used here, with typical S/N mc( i ~ 22 A -1 , mea- 
surements of [Oni] are less reliable than those of [On]. 
In the LINER-like galaxy sample, fewer than one third 
of all galaxies have [Om] EW detections at or above a 
3(7 confidence level whereas more than 75% have > 3cr 
detections in [On]. Because the LINER-like galaxies lie 
on a narrow sequence in EW([On])-EW(Ha) space (see 
Figure 1), [On] can provide an estimate of Ha emission 
which is less dependent upon the Balmer line strengths 
of the stellar template fit used to measure emission lines 
in Yan et al. (2006). Having estimated the Ha EW from 
[Oil] , we can use a theoretically predicted Balmer decre- 
ment to derive the expected emission line strengths of 
H/3, Hy, and R5. 

The existence of a tight roughly linear scaling rela- 
tion between EWQO11]) and EW(Ha) suggests that both 
emission lines are produced by the same mechanism. The 
classification of quiescent vs. LINER-like galaxies was 
based in part on a S/N criterion that both emission 
lines be detected at a 2cr confidence level. It is reason- 
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able to assume that galaxies exist with very faint emis- 
sion that is undetected above the stellar continuum in 
the S/N med ~ 22 A -1 SDSS spectra. These galaxies 
would most likely be intrinsically LINER-likc galaxies 
with emission line equivalent widths approaching zero 
and would be classed by us in the quiescent sample. Us- 
ing a linear least-squares fit of EW(Ha) as a function of 
EWQOii]) and forcing the fit to go through the origin, we 
find that EW(Ha) w 0.197 EWQOii]). This fit is shown 
as the dashed line in Figures la-b. Using this fit, the 
measured EW([On]) for each galaxy is converted into an 
estimated EW(Ha) and then into an estimated Ha flux, 
using the median spectral continuum level around Ha 
in each spectrum. The wavelength range used for com- 
puting the Ha continuum level is the same as Yan et al. 
(2006, Table 3). 

The theoretically predicted Balmcr decrement Ha/H/3 
for LINER- and AGN-like emission line ratios is 3.1 
(Osterbrock & Ferland 2005), with further decrements 
Hy/H/3 = 0.46 and H<5/H/3 = 0.26 for the higher order 
Balmer lines. These values are fairly independent of the 
exact excitation mechanism. However, if there is dust 
within the galaxy in or around the emission line region, 
H/3, Hy and US emission will be dimmed with respect 
to Ha. The Balmer decrement actually measured for 
LINER-like galaxies in the Yan et al. (2006) sample is 
Ha/H/3 = 4.1 ± 0.3, corresponding to reddening with 
r = 0.28 between Ha and H/3, or E(B - V) = 0.47. 

As already noted, the Balmer emission measurements 
in the Yan et al. analysis are highly dependent on 
the choice of stellar template, and Balmer line emis- 
sion is weak in LINER-like objects; H/3 is detected at 
the 3cr level in less than 23% of the galaxies in the 
LINER-like sample. In light of this, we also exam- 
ined the sample of bright galaxies with Bt < 12.5 
mag compiled by Ho, Filippcnko, & Sargent a decade 
ago. (Ho, Filippenko, & Sargent 1995, Ho, Filippcnko, 
& Sargent 1997 and other papers in the series). In 
Ho, Filippenko, & Sargent (1997), they classify emis- 
sion line galaxies based upon their emission line ratios 
(specifically logQOm] A5007/H/3), log([Nn]A6583/Ha), 
and log([Oi]A6300/Ha)) as Seyfert, LINER, Hn, or Tran- 
sition Objects. For all objects in their catalog which are 
classified as LINERs and which satisfy our color cut, the 
median Ha/H/3 Balmer decrement is 3.27, only slightly 
higher than that predicted by theory. This suggests that 
there is very little dust absorption within these systems, 
corresponding to a reddening of r = 0.053 between Ha 
and Rf3, equivalent to E(B - V) « 0.050. The distribu- 
tion of Ha/H/3 in their sample of LINERs is independent 
of galaxy luminosity and emission line strength. How- 
ever, their spectra cover the nuclear region (the central 
~ 200 pc) only, and therefore may not be an appropriate 
counterpart to the SDSS fiber spectra used in our anal- 
ysis, which typically cover 0.35-1.1 galaxy effective radii 
corresponding to several kpc. 

In view of the uncertainty in the value of the proper 
Balmer decrement, the analysis below was repeated for 
both Ha/H/3 = 3.27 (as in Ho, Filippenko, & Sargent 
1997) and Ha/H/3 = 4.1 (as in Yan et al. 2006). To 
obtain an emission infill correction for H7 and H(5, we 
assumed intrinsic Balmer decrements of H7/H/3 = 0.46 
and H<5/H/3 = 0.26, with a reddening law that scales 
as t oc A~ - 7 (following Kauffmann et al. 2003 and 



Cardiel et al. 1998). In §4.1 we will show that the choice 
of Balmer decrement affects the stellar population ages 
measured from the Balmer lines, particularly in H/3. The 
higher order Balmer lines (and H<5 in particular) are less 
affected, and the evidence is strong that while uncertain- 
ties in the Balmer decrement introduce uncertainties into 
our age measurements, they do not qualitatively change 
our main results. 

2.3. Composite Spectra 

Detailed SSP modeling requires high S/N spectra, 
ideally S/N > 100 A" 1 . Individual SDSS spectra in 
this sample typically have S/N mc( i <~ 22 A -1 , mak- 
ing it necessary to stack dozens of similar spectra to 
achieve adequate S/N for SSP modeling. In this pa- 
per, we explore the dependence of stellar population 
age, [Fe/H], and abundance ratios on aperture-corrected 
galaxy velocity dispersion 4 (a) , and the observed LINER- 
like emission line strength as measured by the strength 
of the [On]A3727 emission line. To do this, we divide 
the LINER-like galaxy sample into six bins in cr, with 
cr = 70-120 km s~\ a = 120-145 km s" 1 , a = 145-165 
km s-\ cr = 165-190 km s"\ a = 190-220 km s"\ and 
a = 220-300 km s _1 . Bins in a are chosen to contain 
roughly equal numbers of galaxies. Within each cr-bin, 
the LINER-like galaxies are further divided into bins of 
weak (EWQOii]) < 5 A) and strong (EWQOii]) > 5 A) 
emission strength. Each bin contains over 100 galaxies 
and has a S/N med > 230 A -1 . In total, there are 2141 
galaxies with weak emission and 1775 with strong emis- 
sion. 

Once the LINER-like galaxies have been divided into 
weak and strong subsamples, there are almost five times 
the number of quiescent galaxies as there are either weak 
or strong LINER-like galaxies. For comparison with the 
weak and strong LINER-like galaxies, we construct a 
sample of 2000 quiescent galaxies. Figure 3 shows the 
distribution of EWQOii]) for the 10,284 in the quiescent 
sample as the open histogram along with a gaussian fit 
to the distribution. The distribution peaks <~ 0.3 A be- 
low zero. This may mean that the zeropoint correction 
determined in Yan et al. (2006) is uncertain at this level. 
For LINER-like emission with EWQOii]) ~ 5x EW(Ha) 
and Ha/H/3 ~ 3, this produces AH/3 ~ 0.02 A. An uncer- 
tainty at this level is insignificant in the stellar popula- 
tion modeling process, thus we ignore the possible small 
discrepancy in zeropoint. 

To create a comparison sample of quiescent galax- 
ies, we construct a gaussian distribution centered about 
EWQOii]) = A with a [on] = 1.56 A to match the width 
of the distribution in Figure 3. This distribution is trun- 
cated at ±2(j[oii] in order to exclude outliers which may 
contain low-level undetected emission. From this trun- 
cated gaussian, we randomly select 2000 galaxies for our 
analysis. The EWQOii]) distribution of the final quies- 
cent sample is shown as the shaded histogram in Figure 
3. 

4 Following Bcrnardi ct al. (2003a) and j0rgcnscn, Franx, 
& Kjacrgaard (1995), velocity dispersions are aperture corrected 
to a standard one-eighth of the effective radius using a = 
cr fib( , 'fib/i r o) ' 04 , where is the velocity dispersion measured 
in the fiber spectrum, r is the circular galaxy radius in arc seconds 
and rflb is the spectral fiber radius, 1.5" for SDSS. This correction 
is small, of order ~6%. 
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The 2000 quiescent galaxies are then divided into six 
bins in a as were the weak and strong LINER-like galax- 
ies. The sample therefore includes 18 bins in total: six 
bins in a by three bins in emission line strength, labelled 
as "quiescent", "weak", and "strong". The parameters 
of the various galaxy bins are summarized in Table 1. 

In each individual spectrum, absorption lines are 
broadened by the internal velocity dispersion of the stars 
in the galaxy, which has the effect of reducing the depths 
of the lines in the higher-tr galaxies. To make a consistent 
comparison between separate bins in a, all the spectra 
are smoothed to an effective a = 300 km s^ 1 to match 
the highest-er galaxies in the sample. 5 The spectra are 
coadded within each bin, using an algorithm which av- 
erages the unsmoothed spectra within small sub-bins in 
Ufib before smoothing to a — 300 km s _1 and coadding 
the sub-bins. This minimizes the effect of far outlier 
pixels by averaging out noise spikes in a single pixel be- 
fore smoothing the spectra. The regions around bright 
skylines at 5577A, 6300A, and 6363A (observed frame) 
are masked in the individual galaxies to produce clean 
stacked spectra. In addition to an average spectrum, the 
total S/N at each resolution element is computed to pro- 
duce an error spectrum for each stacked spectrum. 

It is possible, even likely, that the quiescent galaxy 
bins in fact contain galaxies which do have faint emis- 
sion lines but the galaxy spectra are too low S/N for 
the emission lines to be detected at the > 2a level that 
would qualify them for the LINER-like classification. We 
have attempted to remove some of these cases by truncat- 
ing the EW([On]) distribution at ±2<7roii], as described 
above. If a substantial number of weak emission line 
objects remain, these could produce non-negligible emis- 
sion infill in the stacked spectra of quiescent galaxies. A 
check against this is to measure the [On] emission in the 
stacked quiescent galaxy spectra, which have extremely 
high S/N. This is done in the same way as for the in- 
dividual spectra, subtrating off a stellar absorption line 
template before making an emission line measurement, 
using the procedure described in the appendix of Yan et 
al. (2006). The EW([0il]) values measured for the quies- 
cent stacked spectra are -0.01 A, 0.04 A, 0.41 A, -0.27 A, 
0.28 A, and -0.44 A for the a = 70-120 km s"\ 120-145 
km s"\ 145-165 km s" 1 , 165-190 km s" 1 , 190-220 km 
s _1 , and 220-300 km s _1 stacked spectra, respectively. 
The largest of these [On] EWs, 0.41 A, corresponds to 
an Ha EW of 0.082 A and an H/3 EW of 0.027 A, us- 
ing the scaling of equation 2 and the theoretical Balmcr 
decrement. This level of contamination in H/3 is within 
the measurement errors of the H/3 Lick index for these 
high S/N spectra (see §3.1). 

The stacked spectra with 70 km s _1 < a < 120 km s _1 
are shown in Figure 4 as examples of the quality and res- 
olution of the stacked spectra. The black line shows the 
quiescent galaxy bin, the green line shows the galaxies 
with weak LINER-like emission, and the red line shows 
the galaxies with strong LINER-like emission. The weak 
and strong LINER-like bins have been corrected for emis- 

5 Note: Galaxies are binned by aperture-corrected a, as this 
is the best estimate of an intrinsic property of the galaxy. For 
the purpose of measuring absorption line indices, all galaxies are 
smoothed based on their observed to simulate the effect of 
uniformly observing all galaxies with broadening of 300 km s . 



sion infill. The weak and strong bins without the emission 
infill correction are over plotted as dotted lines; the only 
differences appear in the Balmer absorption lines where 
the emission infill correction has been made. Prominant 
emission features are indicated by vertical dashed lines 
and are labelled. Other than the emission features, the 
spectra are remarkably similar, which shows both the 
high S/N of the stacked data and the subtleness of the 
spectral differences that are studied in this paper. 

2.4. Sample Completeness 

The galaxy sample presented here is a magnitude- 
limited sample, as described in §2.1, chosen within a nar- 
row range in redshift (0.06 < z < 0.08). The galaxy bins, 
however, are defined in terms of a. This means that the 
galaxy bins will be incomplete at faint magnitudes, where 
galaxies of a given a may exist but be too faint to make 
it into the sample. Because galaxy magnitude and ve- 
locity dispersion are related through the Faber- Jackson 
relation (Faber & Jackson 1976), the lowest-cr bins will 
be most affected by this incompleteness. If the stellar 
populations of fainter galaxies differ systematically from 
those of brighter galaxies at the same a, the results ob- 
tained from the stacked spectra will be biased at low a. 
There is every reason to expect that the stellar popula- 
tions will indeed differ, since fainter magnitudes at the 
same a likely imply larger mass-to-light ratios. The in- 
complete low-cr bins will be biased toward brighter, and 
therefore possibly younger and/or more metal-poor stel- 
lar populations. 

The top panel of Figure 5 shows a color-magnitude di- 
agram (CMD) of the quiescent and LINER-like galaxy 
samples. Galaxies have been selected by the color cut 
described in §2.1, so only red sequence galaxies are in- 
cluded (compare to Figure lc). The data points are color- 
coded by a as explained in the caption. The dashed 
line at 0A M r = —19.89 shows the completeness limit at 
z = 0.08. The solid line at 01 M r = -19.63 shows the 
magnitude at which the sample is incomplete at the 50% 
level (i.e. missing half the galaxies at that magnitude). 
Thus the dashed line shows where the sample begins to 
miss galaxies and the solid line shows a point at which 
the sample has become severely incomplete. 

In the lower panels of Figure 5, CMDs for each of the 
six ranges in a are shown separately. In each CMD, 
the gray points show the entire sample of quiescent and 
LINER-like galaxies. Overplotted are colored contours 
showing the distribution of galaxies within a given range 
of a. Solid and dashed lines indicate the 50% and 100% 
completeness limits respectively, as in the top panel. Be- 
low each CMD is a histogram of the 01 M r values for the 
galaxies in each a range. 

Two well-known trends are immediately visible: higher 
a galaxies are brighter (Faber- Jackson), and higher a 
galaxies are redder, producing the slope of the red se- 
quence. What is worthy of note is that at fixed a, there 
is little if any slope to the color-magnitude relation. This 
is consistent with the claim of Bernardi et al. (2005) that 
the color-magnitude relation is purely a result of color-c 
and magnitudc-er relations which are more fundamental, 
and Figure 5 gives a straight-forward visual confirmation 
of their result. 

By comparing the distribution of the galaxies in each a 
range (both the 2-dimensional distribution in the CMD 
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as shown in the contour plots, and the 1-dimcnsional his- 
togram in M r ), we can determine the level of incom- 
pleteness in each a range. The lowest-cr range, shown 
in purple, is significantly incomplete; the distribution of 

°- 1 M r has not even reached a peak at the completeness 
limit (dashed line) and the turn-over in the histogram 
may be entirely due to incompleteness. In this a range, 
the stacked spectra will be missing the contribution of 
the faintest galaxies and will thus give an estimate of 
mean stellar population properties which is biased to- 
ward brighter galaxies at fixed a. 

To estimate the incompleteness of each a range, we as- 
sume that the highest a range (a = 220-300 km s _1 ) is 
complete and that its distribution in 01 M r represents 
the true distribution in 01 M r for each range. The his- 
togram of 01 M r for this range (shown in red in the lower 
right corner of Figure 5) shows an extended tail to fainter 
magnitudes. This tail is also visible in the a = 190-200 
km s _1 range and is suggested in the a = 165-190 km s _1 
range, thus lending support to the assumption of similar 

01 M r distributions in all a ranges. In each a range, the 
complete (bright) end of the normalized M r distribu- 
tion is fit to the normalized a = 220-300 km s _1 dis- 
tribution and the discrepancy between the current range 
and the a — 220-300 km s _1 range at the faint end 
is used to estimate the fractional incompleteness of the 
lower a range. We find that the missing galaxies consti- 
tute approximately 51%, 28%, 19%, 11%, and 4% of the 
a = 70-120, 120-145, 145-165, 165-190, and 190-220 
km s _1 ranges, respectively. The estimated fractional 
incompleteness for each a range is shown in Table 1. 

This indicates that the lower a ranges are significantly 
incomplete. If there are systematic differences between 
the stellar populations of bright and faint galaxies at 
fixed cr, the average stellar population parameters de- 
rived from the stacked spectra of decreasing a ranges will 
be increasingly biased toward the bright galaxies at fixed 
a. In appendix A we show that in fact the fainter galaxies 
at fixed a have systematically weaker metal absorption 
lines than the bright galaxies, but that their Balmer ab- 
sorption lines are the same. This suggests that the fainter 
galaxies are more metal poor and somewhat older than 
the brighter galaxies at fixed a. We will correct for the 
effect of this trend in our incomplete data in §3.1. 

3. LICK INDEX LINE STRENGTHS 

3.1. Line Strength Measurements 

In each stacked spectrum, we measure absorption line 
strengths of the Lick indices, as defined in Worthey et 
al. (1994) and Worthey & Ottaviani (1997). These in- 
clude the Balmer lines H/3, Hjf, and HJp, as well as 
numerous "metal lines" : five lines dominated by iron ab- 
sorption (Fe4383, Fe4531, Fe5270, Fe5335, and Fe5406), 
and lines which are sensitive to the abundances of other 
elements (Mgi, Mg 2 , Mg b, CNi, CN 2 , G4300, C 2 4668, 
Ca4227, Ca4455, NaD, TiOi, and Ti0 2 ). The Lick in- 
dices are measured using the IDL code Lick_EW, part of 
the EZ_Ages code package described in Graves & Schi- 
avon (in preparation). Lick_EW uses the error spectra 
described in §2.3 to compute the statistical errors in the 
index measurements due to the finite S/N of the stacked 
spectra. The error calculation is based on the formulae 
of Cardiel et al. (1998) as given in their equations (33) 



and (37). In addition to the Lick indices, the 4000 A 
break has been measured using the narrow definition of 
D„4000 given in Balogh et al. (1999), with errors cal- 
culated from Cardiel et al. (1998) equation (38). Lick 
index line strengths for each stacked spectrum are given 
in Tables 2-4. 

Five Lick indices are excluded from this analysis. 
These include the wide Balmer indices Hja and US a, 
which are designed to measure the broad Balmer wings 
of A stars and are not relevant for the analysis of old 
(> 1 Gyr) stellar populations. Also excluded arc the 
reddest iron lines Fe5709 and Fe5782 because they are 
weaker features than the bluer iron lines and therefore 
represent noisier measurements. Finally, we exclude the 
iron line Fe5015 because it is strongly contaminated by 
[Oiii]A5007 emission in the LINER- like galaxies. 

In order to have a uniform sample, the stacked spectra 
are all smoothed to a = 300 km s _1 to match the a of 
the most massive galaxies (see §2.3). The combined reso- 
lution of g == 300 km s _1 and the native SDSS resolution 
is at lower resolution than the Lick/IDS instrument res- 
olution at which the Lick indices are defined. The index 
measurements must therefore be corrected back to effec- 
tive Lick/IDS resolution. This correction is performed 
by Lick_EW using velocity dispersion corrections deter- 
mined from smoothed SSP spectra from Schiavon (2007, 
Table A2a). The indices reported here are therefore ef- 
fectively at Lick/IDS resolution. However, no zeropoint 
correction is performed to place them on the original 
Lick/IDS system. This is for two reasons. Firstly, the 
Lick standard stars are not observed by the SDSS — they 
are too bright — thus the data needed to make this cor- 
rection do not exist. More importantly, the SSP mod- 
els which we use to interpret the index measurements 
in §4.1 are defined using flux-calibrated spectra and are 
themselves not zeropoint-shifted to the Lick/IDS system. 
Because the SDSS spectra are also flux-calibrated, any 
zeropoint offset between the indices measured in SDSS 
spectra and the models should be small (see discussion 
in Schiavon 2007). 

Figure 6 shows the indices measured in the stacked 
spectra as a function of the mean a in each bin. Colors 
correspond to the emission line properties of the galax- 
ies as follows: quiescent galaxy bins are shown in black, 
those with weak LINER-like emission in green, and those 
with strong LINER-like emission in red. In each of the 
panels a-d, the left column shows the Lick index line 
strengths measured in the stacked spectra. The error 
bars include statistical errors due to finite S/N only — 
no systematic errors are included. The lines show least 
squares fits of index strength onto a for each bin in emis- 
sion line strength. 

In §2.4, we demonstrated that the lower a bins are 
significantly incomplete at faint magnitudes. A correc- 
tion for this effect, described in detail in Appendix A, is 
applied to the measured Lick index line strengths; these 
corrected line strengths are shown in the right columns of 
panels a-d in Figure 6. For the moment, we note merely 
that the correction is a small effect in all indices and has 
no qualitative effect on the following discussion. 

No attempt has been made to correct the line strengths 
for variations in aperture. The range of rcdshifts in this 
sample spans a factor of 30% in angular diameter dis- 
tance, but because the redshift distributions of quies- 
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cent, weak, and strong LINER-like galaxies are similar, 
the differences in mean physical aperture radius between 
the samples under comparison are negligible. 

3.2. Line Strength Trends with a and Emission 
Properties 

In SSP models with ages over 1 Gyr, Balmer line 
strengths decrease with increasing SSP age. In the 
same models, the 4000 A break (D„4000), the CH band 
(G4300), and all of the metal line strengths increase 
with increasing SSP age (Worthey et al. 1994, Worthcy 
& Ottaviani 1997, Bruzual & Chariot 2003, Thomas, 
Maraston, & Bender 2003, Schiavon 2007). However, line 
strengths are also affected by the metallicity of the SSP. 
In the same models, increasing the metallicity changes 
the line strengths in the same direction as increasing the 
SSP age: Balmer line strengths decrease while D„4000, 
G4300, and metal line strengths increase with increasing 
SSP metallicity. These effects are summarized in Table 
5 for those unfamiliar with SSP model results. 

Figure 6a shows the Balmer lines H/3, H^f, and HSp, 
as well as D n 4000 and G4300. All of these indices 
vary systematically with a such that, as a increases, 
the Balmer lines get weaker and D„4000 and G4300 get 
stronger. Panels b-d in Figure 6 show line strengths 
for the metal lines. These indices show the same trends 
as D n 4000, and G4300: composite galaxy spectra with 
higher a have stronger metal absorption lines than those 
with low a. The stellar populations of the galaxies pre- 
sented here vary systematically with a, as has been seen 
previously by many authors. 

It is also clear from Figure 6 that, at fixed a, galax- 
ies with strong LINER-like emission (red points) typi- 
cally have stronger Balmer line absorption and weaker 
D„4000, G4300, and metal lines than galaxies with weak 
LINER-like emission (green) , which in turn have stronger 
Balmer lines and weaker D„4000, G4300, and metal lines 
than the quiescent galaxies (black). This suggests that, 
in addition to a relationship between the stellar popula- 
tions and a, stellar populations also vary with the emis- 
sion line strength of the galaxy. This trend exists not 
only in the Balmer lines, which are sensitive to errors in 
the emission infill correction made in §2.2, but also in 
D„4000, G4300, and the metal lines which are not af- 
fected by contamination from either the emission or the 
chosen emission infill correction. 

The sodium line NaD is the one exception to the trend 
that metal lines weaken as LINER-like emission strength 
increases; NaD absorption strengths increase with in- 
creasing LINER-like emission and typically vary by sev- 
eral tenths of an A between quiescent and strong emission 
galaxies at fixed cr. NaD is a resonance line and is known 
to be significantly affected by interstellar absorption. In 
their sample of Galactic globular clusters, Gorgas et al. 
(1993) assume that up to 0.8 A of the observed NaD line 
strengths is due to interstellar rather than stellar absorp- 
tion. If the intrinsic stellar NaD line strength in galaxies 
actually follows the behavior of all the other metal lines 
and decreases with increasing emission, then the observed 
increase of NaD absorption with emission suggests that 
LINER-like galaxies contain more cold interstellar ma- 
terial than their quiescent counterparts, which serves to 
enhance the observed NaD absorption. This is discussed 
in greater detail in §5.1. 



The trends visible in Figure 6 (with the exception of 
NaD, as discussed above) are consistent with the hypoth- 
esis that stars in galaxies with higher a are typically 
older than those in lower-cr galaxies, and that at fixed 
a galaxies with LINER-like emission have younger stars 
on average than their quiescent counterparts. However, 
as shown in Table 5, increasing SSP age changes all in- 
dex strengths in the same direction as increasing SSP 
metallicity This means that the observed trends with 
a and emission strength may be trends in age, trends 
in metallicity, or some combination of the two. For- 
tunately, although the response of the various lines to 
changes in age is qualitatively similar to the response to 
changes in metallicity, there are quantitative differences: 
the Balmer lines are more sensitive to changes in age and 
less sensitive to changes in metallicity than are the metal 
lines. By comparing SSP model predictions to pairs of 
measured line strengths that differ quantitatively in their 
sensitivity to age and metallicity, we can break the "age- 
mctallicity degeneracy." This is done in §4.1. 

3.3. Effect of the Infill Correction on Balmer Line 
Strengths 

Figure 6 shows that galaxies with LINER-like emission 
typically have stronger Balmer absorption lines than qui- 
escent galaxies at the same a. However, the strength of 
this effect depends critically upon the assumptions made 
in §2.2 in correcting for emission infill in these galaxies. 
It is instructive to examine the total effect of the emis- 
sion infill corrections to understand their significance for 
Balmer absorption line strengths. 

Figure 7a-c plots the Balmer lines as measured in the 
composite spectra without any emission infill correction. 
As in Figure 6, black, green, and red data points show 
quiescent, weak, and strong LINER-like emission galaxy 
bins, respectively. The emission infill is a large effect 
in H/3, showing significant infill from emission in both 
weak and strong LINER-like galaxies. However, the infill 
in H"fF is substantially less, and is least significant in 
RSf due to the Balmer decrement. Even without any 
emission infill correction at all, H5f still shows a trend 
of stronger Balmer absorption in stronger LINER-like 
galaxies. Moreover, making no infill correction at all is 
obviously wrong, as can be seen by the substantial infill 
in the uncorrected H/3. Any infill correction will only 
make the trend in H5f stronger. 

Figures 7d-f and 7g-i show the Balmer lines measured 
with the two different assumptions for the Balmer decre- 
ment that were described in §2.2: Ha/H/3 = 4.1, as mea- 
sured by Yan et al. (2006), and Ha/H/3 = 3.27, as seen 
in the Ho, Filippenko, & Sargent (1997) red sequence 
LINERs. The choice of Balmer decrement makes some 
difference in the H/3 measurement (as can be seen com- 
paring Figure 7f and 7i), but it is a small effect for Rjf 
(Figure 7e and 7h) and is practically negligible for H<5f 
(Figure 7d and 7g). 

The infill correction Ha/H/3 = 4.1 gives H/3 absorption 
line strength values that cross at high cr; at a < 190 km 
s _1 galaxies with LINER-like emission have stronger H/3 
absorption than their quiescent counterparts, but above 
cr = 220 km s _1 , the H/3 absorption in the LINER-like 
galaxies is weaker than in the quiescent galaxies. In §4.1 
we will show that this results in age trends that reverse 
for galaxies with cr < 190 (where LINER-like galaxies are 
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younger than the quiescent galaxies) versus those with 
a > 220 (where LINER-like galaxies are older than the 
quiescent galaxies). It is difficult to find a physical jus- 
tification for this reversal of the trend in H/3 (and there- 
fore age) with emission strength, and furthermore this 
reversal is not corroborated by Hjf or RSp; the sim- 
pler explanation is that assuming a Balmer decrement of 
Ha/H/3 = 4.1 underestimates the appropriate emission 
infill correction. Figure 7f shows that the reversal de- 
scribed here is hardly a significant effect but would be 
larger for even larger values of the Balmer decrement; 
therefore this argument suggests that Ha/H/3 = 4.1 is 
an upper limit on the Balmer decrement. 

Using Ha/H/3 = 3.27 gives consistent trends with 
LINER-like emission strength versus a for all Balmer 
lines, with H/3 line strengths that converge to a single 
value for all emission strengths at high a. If the in- 
trinsic Balmer decrement for these LINER-like galaxies 
is similar to that of AGN (which have Ha/H/3 = 3.1, 
Osterbrock & Ferland 2005), then this is very close to 
the theoretical (unreddened) lower limit on the Balmer 
decrement. It is therefore likely that these two values 
enclose the true mean value, and 3.27 < Ha/H/3 < 4.1 
for these systems. 

The assumption that LINER-like emission is similar to 
AGN emission is not necessarily a legitimate one. Emis- 
sion from stellar- ionized Hn regions has Ha/H/3 = 2.86 
in the absence of dust absorption (Osterbrock & Ferland 
2005). For the sake of completeness, Figure 7 also shows 
the Balmer lines with an infill correction assuming this 
theoretical Hn region value. The strength of H/3 changes 
noticeably, but the effect on H<5p is within the measure- 
ment errors. In summary, Ha/H/3 = 3.27 appears to be 
a reasonable value for the mean reddened Balmer decre- 
ment and it is almost certain that 2.86 < Ha/H/3 < 4.1. 
From Figure 7, it is clear that the exact value of the 
mean Ha/H/3 has some effect on H/3 line strengths but 
does not qualitatively changing the trends described in 
§3.2 above, and that H^f and H<5f are only marginally 
affected by the exact value of the decrement. 

4. STELLAR POPULATION MODELING: 
PROCEDURES AND CAVEATS 

4.1. SSP Models 

Stellar population models make it possible to convert 
measured index strengths into fundamental stellar popu- 
lation parameters (i.e., age, [Fe/H], and elemental abun- 
dance ratios). In this paper, we have chosen to use 
the models of Schiavon (2006; hereafter S06), which are 
based on fitting functions. Although the models produce 
SSPs rather than more complicated star formation his- 
tories, they include the effects of multiple individual ele- 
mental abundances in addition to mean age and [Fe/H]. 

The S06 models calculate predicted Lick index values 
for a model SSP with a specified abundance pattern. 
Non-solar abundance patterns change the SSP model 
through two different channels. The first is the effect on 
the isochrone of a non-solar abundance pattern, which 
changes the temperature and luminosity of the main se- 
quence turn-off and red giant branch. The S06 models 
incorporate this effect (although in a limited way) by 
providing the option to construct the SSP model based 
on one of two isochrones: the solar-scaled isochrone of 



Girardi et al. (2000) or an a-enhanced isochrone from 
Salasnich et al. (2000) with average [a/Fe] = +0.42. 6 
The second way in which non-solar abundances enter the 
SSP models is through their effect on the various line 
opacities of stellar atmospheres. In modeling line opaci- 
ties, the S06 models take into account the effects of the 
elements C, N, O, Mg, Ca, Na, Si, Cr, and Ti. Stellar 
line opacities due to various elements and molecules may 
lie within the index passband of a Lick index, enhancing 
the strength of the measured index. They may also lie in 
the blue and red continuum regions on either side of the 
index against with the index absorption strength is mea- 
sured. In this latter case, increased absorption weakens 
the measured Lick index strength. The Korn, Maraston, 
& Thomas (2005) index response functions are used to 
calculate the effect on the model Lick indices due to each 
of the elements listed above. 

It should be emphasized that the S06 models are com- 
puted at fixed [Fe/H], rather than at fixed total metallic- 
ity ([Z/H]), as is the case in many other stellar popula- 
tion models. This has several benefits. The first is one of 
simplicity and conceptual transparency. For models com- 
puted at fixed [Z/H], changing the a-enhancement of the 
model also changes [Fe/H]. Total metallicity is in fact 
dominated by oxygen (an a-elemcnt), so that changing 
[a/Fe] a small amount results in a large change in [Fe/H] 
in order to keep [Z/H] constant. It is therefore concep- 
tually more complicated to keep track of how enhancing 
[a/Fe] affects stellar populations because one must dis- 
entangle the effects of increasing a-elements from the 
effects of decreasing [Fe/H]. One case where this is par- 
ticularly important is that of the higher order Balmer 
lines, which are strongly affected by Fe I absorption 
lines in their continuum (Thomas, Maraston, & Korn 
2004). In fixed [Z/H] models, changing [a/Fe] dramat- 
ically changes [Fe/H], which in turn has a strong effect 
on Hjf and US p. Because the S06 models are computed 
at fixed [Fe/H], the model predictions for Hyp and KSp 
are relatively insensitive to [a/Fe]. This is discussed in 
detail in §5.2. 

A second benefit of computing models at fixed [Fe/H] 
rather than fixed [Z/H] is that oxygen, which dominates 
total metallicity, is notoriously difficult to measure in 
the integrated light of stellar populations. While [O/H] 
remains essentially unmeasured, it is safer not to claim 
knowledge of total [Z/H]. This will be discussed in more 
detail in §4.3 and §5.3. A related point is that [Fe/H] is 
a parameter that is more directly related to measurable 
quantities (i.e., iron absorption lines such as Fe5270 or 
Fe5335) than is [Z/H], which depends on the abundances 
of many different elements. Some of these elements are 
not measured by any stellar absorption line, including 
the dominant element oxygen. 

The S06 models take as input a set of abundance ratios 
and calculate Lick index values for the given abundances 
at a range of ages and [Fe/H]. The model produces grids 
of Lick index values for ages between 1.2 and 17.7 Gyr, 
and —1.3 <[Fc/H]< +0.2 for the solar-scale isochrone 
and —0.8 <[Fc/H]< +0.3 for the a-enhanced isochrone. 

6 The Salasnich et al. (2000) a-enhanced isochrone was com- 
puted using the following abundances: [O/Fc] = +0.50, [Ne/Fe] = 
+0.29, [Mg/Fe] = +0.40, [Si/Fe] = +0.30, [S/Fc] = +0.33, [Ca/Fe] 
= +0.50, [Ti/Fc] = +0.63, [Ni/Fe] = +0.02. All other abundances 
are solar. 
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The model grids can be used with index-index plots in 
the style pioneered by Worthey et al. (1994) to determine 
the age and [Fe/H] of an observed stellar population. Be- 
cause the various Lick indices have different sensitivities 
to each elemental abundance, a best-fitting set of abun- 
dances can be found by tweaking the input abundances 
in the S06 models until all of the index-index plots give 
consistent values for the age and [Fe/H] of the stellar 
population in question. 

The S06 models are limited to a subset of the Lick 
indices due to the restricted wavelength coverage of the 
stellar spectra used to construct the fitting functions. 
Only the Balmer indices (H/3, along with both the broad 
and narrow H7 and R5 indices), Fe4383, Fe5015, Fe5270, 
Fe5335, Mg 2 , Mg 6, G4300, C 2 4668, CNi, CN 2 , and 
Ca4227 are modeled. This set of indices makes it possi- 
ble to constrain the abundances [Mg/Fe], [C/Fe], [N/Fe], 
and [Ca/Fe], but [O/Fe], [Na/Fe], [Si/Fe], [Cr/Fe], and 
[Ti/Fe] are relatively unconstrained. In our analysis, 
we predominantly use the solar-scaled isochrone and 
set [O/Fe] to match the isochrone (i.e. [O/Fe] = 0.0). 
We also show results computed using the a-enhanced 
isochrone and with [O/Fe]=+0.3 for comparison. Be- 
cause Na, Si, and Ti are a-elements, we force their abun- 
dances to follow [Mg/Fe], while the iron-peak element Cr 
tracks Fc ([Cr/Fe] = 0.0). 

4.2. Fitting for Ages and Abundances 

Fitting for the best set of input abundance ratios is 
done using the IDL code package EZ_Ages (Graves & 
Schiavon, in preparation). Briefly, EZ_Ages implements 
the method developed by Schiavon (2007) to determine 
abundances and ages of stellar populations. It uses an 
index-index plot of H/3 and (Fe) (defined by (Fe) = 
0.5(Fe5270 + Fc5335)) to measure a "fiducial" age and 
[Fe/H]. These indices are chosen because they are rela- 
tively insensitive to other, non-solar abundance ratios. 
EZ_Ages then uses an index-index plot of H/3 and Mg 
b to constrain [Mg/Fe]. This is done by varying [Mg/Fe] 
as input to the S06 model until the age and [Fe/H] mea- 
sured in the H/3-Mg b plot match the fiducial values from 
the H/3-(Fe) plot. In similar fashion, C 2 4668 is used to 
determine [C/Fe], CN 2 is used to determine [N/Fc] once 
[C/Fe] is set, and Ca4227 is used to get [Ca/Fe]. The 
entire process is then iterated in order to produce a self- 
consistent set of abundances. EZ_Ages then measures 
separately an age from Hjf an d an age from H<5p using 
the best-fitting abundances. The result is a measurement 
of [Fe/H], [Mg/Fe], [C/Fe], [N/Fe], and [Ca/Fe], as well 
as ages measured from each of the three Balmer lines, for 
each set of measured Lick indices. 

This abundance fitting must be done in a proscribed 
order. This takes into account the complicated effects of 
the various element abundances on the Lick indices used 
in the fitting procedure. The goal is to only introduce 
one new abundance ratio at a time. The Mg b index is 
sensitive to age, [Fe/H], and [Mg/Fe], but is relatively 
insensitive to all other abundance ratios (Korn, Maras- 
ton, & Thomas 2005, Serven, Worthey, & Briley 2005). 
Therefore, once the fiducial age and [Fe/H] have been 
determined, Mg b can be used to fit for [Mg/Fe]. The 
abundance ratios [C/Fe], [N/Fe], and [Ca/Fe] are more 
difficult to isolate. The [N/Fe] abundance ratio is diffi- 
cult because its effect is seen most strongly in the CNi 



and CN 2 indices where [C/Fe] and [O/Fe] are also impor- 
tant factors; the Ca4227 index is highly sensitive to CN, 
as well as being sensitive to [Ca/Fe] because of a strong 
CN band in the blue continuum region around the in- 
dex (described in detail in Prochaska, Rose, & Schiavon 
2005); and [C/Fe] can best be isolated in the C 2 4668 in- 
dex, however this index is also sensitive to [O /Fe] (Korn, 
Maraston, & Thomas 2005, Serven, Worthey, & Briley 
2005). As mentioned above, oxygen is very difficult to 
measure in unresolved stellar populations. 

In the EZ_Ages abudancc fitting process, [C/Fe] is 
set first by fitting the C 2 4668 index, because C 2 4668 de- 
pends only on age, [Fe/H], [C/Fe], and [O/Fe]. ThenCN 2 
can be used to fit [N/Fe] because [C/Fe] is already de- 
termined. Only once [C/Fe] and [N/Fe] are determined 
can Ca4227 be used to fit for [Ca/Fe], in order to ac- 
count for the strong dependence of Ca4227 on CN. The 
fact that [O/Fe] remains unknown contributes a signif- 
icant uncertainty to all of these abundance determina- 
tions, but is unavoidable. When presenting abundance 
measurements for our data in §5.3, we will illustrate the 
effect of changing [O/Fe] by 0.3 dex. In practice, chang- 
ing [O/Fe] by this amount appears to have no effect on 
[Fe/H] or [Mg/Fe], some effect on [C/Fe], a substantial 
effect on [N/Fe], and almost no effect on [Ca/Fe]. 

In order to illustrate the effect of elemental abundances 
on different line indices and our abundance-fitting pro- 
cess, we show in Figures 8 and 9 a set of index-index 
plots for all of the Lick indices included in the S06 model. 
In each plot, solid lines connect models at fixed [Fe/H], 
while the dotted lines connect models of constant age. 
The three data points are from the stacked spectra of 
quiescent galaxies with a — 70-120 km s _1 , a = 145- 
165 km s" 1 , and a = 220-300 km s" 1 plotted as the 
triangle, star, and square, respectively. The index mea- 
surement errors are smaller than the size of the plotting 
symbols. 

The model shown here is computed with the EZ_Ages 
output best-fitting abundances for the a — 70-120 km 
s- 1 stacked spectrum: [Mg/Fe] = 0.120, [C/Fe] = 0.100, 
[N/Fe] = -0.002, and [Ca/Fe] = 0.010, calculated with 
a solar-scale isochrone, with solar [O/Fe], with Na, Si, 
and Ti abundances set to equal [Mg/Fe], and with so- 
lar [Cr/Fe]. Figure 8 shows the Lick indices which are 
used by EZ_Ages in the abundance fitting. Note that 
the a = 70-120 km s -1 data point (triangle) which was 
used in the abundance fitting falls in the same place with 
respect to the grid lines in each index- index plot. This is 
the essence of the abundance fitting process and signifies 
that a mutually self-consistent set of stellar population 
parameters and abundance ratios has been found. 

A check on the robustness of the abundance results is 
to compare these index-index plots with those in Figure 
9, which were not used in the abundance fitting process. 
The location of the a = 70-120 km s -1 data point in 
each grid in Figure 9 is consistent with its position in 
the grids of Figure 8, except for the G4300 index which 
is too strong in the best-fitting model and H7^ which is 
too weak in the best-fitting model. A similar discrepancy 
in G4300 was noted in Schiavon (2007) in comparisons 
with Galactic clusters M67 and NGC 6528. Sanchez- 
Blazquez et al. (2006) also show anomalous behavior in 
the G4300 index with respect to C 2 4668 in their sample 
of early- type galaxies. The G4300 index in globular clus- 
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ters appears to plateau for iron abundances with [Fc/H] 
> —0.7; above this [Fe/H] the G band feature contin- 
ues to increase with metallicity but absorption lines in 
the continuum against which the G4300 line strength is 
measured also increase in strength in such a way as to 
maintain relatively constant values of G4300 over a range 
of metallicity. The C24668 index is in a region of the 
spectrum where it is possible to define continua which 
are less affected by absorption. We thus consider C24668 
to be a more reliable indicator of carbon enhancement 
than G4300. As for Hjp, we will show in §5.2 that H-jf 
and H(5p generally give younger age estimates than H/3 
based on a SSP model and discuss the implications of 
this effect. 

Unlike the a = 70-120 km s _1 data, the a = 145- 165 
km s _1 and a = 220-300 km s _1 data points (star and 
square) do not fall in the same location on each grid. The 
Mg, CN, and C24668 lines are particularly discrepant. 
This is because the grids were constructed to match the 
a = 70-120 km s _1 data. In order to produce model 
grids that match these higher a data points, [Mg/Fe], 
[C/Fe], and [N/Fe] must be increased (which shifts the 
model grids for the Mg, CN, and C24668 indices in Fig- 
ures 8 and 9 to the right, toward the data). This suggests 
that higher a galaxies have larger [Mg/Fe], [C/Fe], and 
[N/Fe]. In §5.3 we show that this is indeed in the case. 

For comparison, Figures 10 and 11 show index- index 
plots computed with a different set of abundance ratios. 
Here we have used the best-fitting abundances for the 
220 < a < 300 km s" 1 data, which have [Mg/Fe] = 
0.260, [C/Fe] = 0.270, [N/Fc] = 0.220, and [Ca/Fe] = 
0.040. As before, the other a-elements (Na, Si, and Ti) 
follow Mg, while O and Cr have solar abundances. In 
Figure 10, the 220 < a < 300 km s _1 data (squares) 
give the same age and [Fe/H] values in each index-index 
plot. Because this model has enhanced [Mg/Fe], [C/Fe], 
and [N/Fe], the plots of H/3 against Mg 6, C 2 4668, and 
CN2 are now consistent with the H/3-(Fe) plot for the 
220 < a < 300 km s _1 data. In contrast, the lower a 
data points in these plots lie at values of [Fe/H] that are 
too low to match those from the H/?-(Fe) plot. Figure 
11 shows the indices that are not used in the abundance 
fitting. As before, the plots of H/3 against Fe4383, Mg2, 
and CNi are consistent for the 220 < a < 300 km s _1 
data, while the G4300 index produces ages that are too 
young. Here the ages measured from Hjf and RSf are 
somewhat younger than those measured from H/3. 

An important caveat applies to the ages and abun- 
dances presented in the following sections. The S06 
model fitting functions are constructed from flux- 
calibrated stellar spectra with no zeropoint correction to 
bring them onto the original Lick/IDS system. The SDSS 
spectra used in this analysis are also flux-calibrated, so 
zeropoint offsets between the model index system and 
the data should be small. However, as with most stel- 
lar population models based upon index measurements, 
relative values of stellar population parameters are more 
reliable than absolute measurements due to the zeropoint 
uncertainties. 

4.3. Oxygen Abundances in Old Stellar Populations 

It is unfortunate that it has not been possible thus 
far to measure oxygen abundances in unresolved stellar 
populations. Not only is oxygen the single most abun- 



dant metal in the Universe, but main sequence turn-off 
temperatures are substantially affected by oxygen abun- 
dances (Vandenberg & Bell 2001). An inability to de- 
termine the oxygen abundance therefore means that the 
total metallicity is unknown and the true isochrone shape 
is not well-known, raising uncertainties in stellar popu- 
lation analysis. 

Theoretical models of nucleosynthesis predict that Mg 
and O both originate in Type II SNe (e.g. Wheeler, Sne- 
den, & Truran 1989, Woosley & Weaver 1995) and should 
track one another closely, but there is some evidence that 
this is not always the case. In a survey of 27 Milky Way 
bulge giants, Fulbright, McWilliam, & Rich (2006) find 
that, although Mg appears to be enhanced at all [Fe/H], 
[O/Fe] declines with increasing [Fe/H] and is solar or 
mildly sub-solar at [Fe/H] > 0.0. This indicates that, 
at least in some stellar populations, the assumption that 
O varies in lockstep with Mg is not valid. In addition, 
there is growing evidence from X-ray observations that 
the hot gas in elliptical galaxies has subsolar values of 
[O/Fe] (Humphrey & Buote 2006, Humphrey, Buote, & 
Canizares 2004 and references therein) and also subsolar 
[O/Mg] (Humphrey & Buote 2006), further undermining 
the assumption that Mg and O vary together. 

In light of these uncertainties, we treat [O/Fe] as un- 
determined. The S06 SSP models are computed for fixed 
[Fe/H] rather than fixed [Z/H] so the uncertainty in to- 
tal metallicity should not have a large effect on the re- 
sults of the stellar population modeling. In the follow- 
ing sections, we present results based on both the so- 
lar abundance and a-enhanced isochrones available with 
EZ_Ages, and at each stage of the analysis illustrate the 
effect of changing the oxygen abundance by exploring a 
range of values for [O /Fe] . 

5. STELLAR POPULATION MODELING RESULTS: 
AGES AND ABUNDANCES 

5.1. SSP Ages as a Function of cr and Emission Line 

Strength 

The EZ_Ages code was used to estimate mean 
luminosity-weighted ages and abundances from Lick in- 
dex measurements for each of the stacked spectra. Us- 
ing the best-fitting abundance models, we obtained ages 
for each stacked spectrum, measured separately from the 
H/?-(Fe), H 7F -(Fe), and H<5p-(Fe) grids. We shall here- 
after refer to these separate age measurements as "H/3 
ages", "H7F ages", and u H5f ages", with the under- 
standing that these indicate the ages as measured from 
index-index plots using the named Balmer line. Figure 
12 shows these age measurements as a function of the 
mean a in each galaxy bin and emission line strength, 
for each of the stacked spectra. The left panel shows the 
measured ages when a Balmer decrement of Ha/H/3 = 
3.27 is assumed, while the center panel assumes Ha/H/3 
= 4.1. In both cases, the EZ_Ages fitting has been com- 
puted with a solar-scaled abundance isochrone. The right 
panel shows the effect of using the a-enhanced isochrone 
and [O/Fc] = +0.3 in the EZ_Ages fitting process. The 
H/3 ages, H7F ages, and H<5p ages are plotted separately 
in each panel, from top to bottom. The error bars show 
the uncertainty in the ages due to measurement errors in 
the Balmer lines and (Fe) only. The solid lines are linear 
least squares fits of age onto cr. 
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The most striking result in this figure is that the galax- 
ies with strong LINER-likc emission are systematically 
younger than their quiescent counterparts at the same 
a, with weak LINER-like galaxies falling in between. 
This trend exists in all of the age measurements, regard- 
less of which Balmer line is used, and for either Balmcr 
decrement. As expected from the line strength variations 
shown in Figure 7, the H/3 ages are most affected by the 
chosen decrement, while the H5p ages are essentially un- 
affected. In §3.3, we argued that the two values of the 
Balmer decrement used here likely represent lower and 
upper boundaries for reasonable estimates of Ha/H/3. 
Thus the panels of Figure 12 show the range of possi- 
ble SSP age differences between quiescent and LINER- 
like galaxies. The use of an a-enhanced isochrone with 
[O/Fe] = +0.3 (right panel) results in H/3 ages that are 
consistently younger than those measured using a solar 
abundance isochrone but does not change the trend in 
age with emission strength. 

The difference in measured age between quiescent and 
strongly emitting galaxies is substantial — ranging from 
around 3.5 Gyr in the lowest a bin to 2 Gyr in the highest 
a bin with Hoj/H/3 = 3.27 and a solar-scaled isochrone, 
as measured in all Balmer lines. Galaxies with LINER- 
like emission contain stars which are on average 10-40% 
younger than those in quiescent galaxies of the same a 
(using H/3 ages). 

A second effect that is apparent in Figure 12 is that 
galaxies with high a are on average older than those with 
low a. This trend is apparent in ages measured by all 
Balmer lines, independent of the chosen Balmer decre- 
ment, and for quiescent as well as emission line galaxies 
(with the exception of H<5p ages in quiescent galaxies 
which are ~ 7 Gyr, regardless of a). This suggests that 
more massive galaxies formed a larger percentage of their 
stars in the distant past than did less massive galaxies, 
which seems to be consistent with the so-called "down- 
sizing" phenomenon (Cowie et al. 1996, see also Bernardi 
et al. 2003b, Nelan et al. 2005, Thomas et al. 2005, Trcu 
et al. 2005, and many others). 

Interestingly, the Schiavon (2007) analysis of SDSS 
stacked spectra from Eisenstein et al. (2003) shows no 
increase in H/3 age with galaxy absolute magnitude but 
a definite increase in US age, using the same models as 
in this work. That sample is selected predominantly by 
morphology and color, with no selection applied based on 
emission line properties. Because of this, the Eisenstein 
et al. (2003) sample likely contains not only quiescent and 
LINER-like galaxies, but a small contaminating fraction 
of galaxies with Seyfert and Hil region emission as well. 
The emission infill correction applied by Schiavon (2007) 
is inaccurate for these galaxies, and appears to cause 
the discrepancy in the trends seen in Schiavon (2007) 
with the results presented here. Appendix B presents 
a detailed comparison between Schiavon (2007) and this 
work. 

The stellar population modeling described above shows 
that the LINER-like galaxies have younger ages than 
their quiescent counterparts at fixed <r; it is natural to 
assume that they should therefore have slightly bluer col- 
ors. In fact, the opposite is true. As shown in Figure 13, 
the LINER-like galaxies typically have slightly redder col- 
ors than the quiescent galaxies at the same a. The effect 
is small, ranging up to 0.037 mag in 01 (g-r) but is consis- 



tent across all galaxy bins. If the LINER-like galaxies are 
indeed younger than the quiescent galaxies at the same 
cr, they must also be subject to a larger degree of intrin- 
sic interstellar reddening, on average. This hypothesis is 
corroborated by the observation that LINER-like galax- 
ies appear to have stronger NaD absorption than their 
quiescent counterparts, suggesting that they are subject 
to a greater amount of internal interstellar absorption 
(see the discussion in §3.2). 

More quantitatively, the younger ages measured for the 
LINER-like galaxies in the spectral line fitting (Figure 
12) should make the LINER-like galaxies ~ 0.06 mag 
bluer in B — V , based on the stellar population mod- 
els of Bruzual & Chariot (2003). Instead, they are ob- 
served to be ~ 0.04 mag redder in B — V J This implies 
that the typical LINER-like galaxy is subject to inter- 
nal reddening of E{B — V) w 0.10 beyond the (typi- 
cally small) internal reddening of quiescent red sequence 
galaxies. This modest amount of reddening is reason- 
ably consistent with the value E(B — V) = 0.05 that was 
used in our preferred computation of the Balmer decre- 
ment in §2.2, and is substantially smaller than the value 
E(B — V) = 0.47 that gives an upper boundary to rea- 
sonable values of Ha/H/3 above. 

In the Milky Way, E(B - V) = 0.10 in the Galactic 
plane roughly corresponds to interstellar NaD absorption 
of ~ 0.35 A (Gorgas et al. 1993). This is generally con- 
sistent with the excess of NaD absorption in the LINER- 
like galaxies needed to bring the NaD stellar absorption 
strengths into line with the rest of the metal lines shown 
in Figure 6. The relation between E(B — V) and NaD 
interstellar absorption observed in the Milky Way varies 
substantially outside the Galactic plane, making a calcu- 
lation of E(B — V) from observed NaD absorption nec- 
essarily inaccurate. However, it is clear that there is 
at least reasonable agreement between the internal red- 
dening suggested by the median galaxy colors and the 
observed interstellar NaD absorption in the LINER-like 
galaxies. Thus, LINER-like galaxies typically suffer from 
a greater degree of internal reddening than their quies- 
cent counterparts at fixed a but this excess reddening is 
small, around E(B — V) = 0.1 on average. 

5.2. Ages Measured by Hjf and H5p 

A third effect that is visible in Figure 12 is that ages 
measured with different Balmer lines vary systematically, 
such that bluer Balmer lines produce younger age mea- 
surements. This is true for the quiescent galaxies as well 
as the LINER-like galaxies, so it cannot be an artifact of 
the emission infill correction. 

The Balmer lines Yi-yp and H5f have been shown to 
be sensitive to non-solar abundance ratios in models 
computed at fixed [Z/H]. The effect of non-solar abun- 
dance ratios is strongest at high metallicity (Thomas, 
Maraston, & Korn 2004), and is thus particularly rele- 
vant for early type galaxies, whose metallicities are typ- 
ically nearly solar or super-solar, unlike globular clus- 
ters, which are typically a-enhanced but are compari- 
tively metal-poor. Thomas, Maraston, & Korn (2004) 

7 The 0,1 (g— r) color has been converted to rest frame B—V using 
Mike Blanton's kcorrect code, v. 3_2 (Blanton et al. 2003a) along 
with the empirical transformations from the SDSS ugriz magni- 
tude system to Johnson B and V bands given by Jordi, Grebel, & 
Ammon (2006). 
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have demonstrated that, in models computed at fixed 
[Z/H], including the effects of a-enhanccment raises the 
age estimates from R-yp and H£p, and could thus in prin- 
cipal solve a discrepancy in which the bluer Balmer lines 
are producing younger age estimates than H/3. 

However, the sensitivity of Hjf and H5f is primarily 
a sensitivity to [Fe/H], rather than [a/Fe], and is an ex- 
ample of the complication of computing models at fixed 
[Z/H]. The abundance ratio dependence oiHjF and HS F 
is primarily due to the large number of iron absorption 
lines in the regions surrounding the Hjf and H5f indices 
(Thomas, Maraston, & Korn 2004). For models at fixed 
[Z/H], choosing a model with too low a value of [a/Fc] 
(i.e., a solar-scaled model when the galaxy is a-enhanced) 
means that the model has too high a value of [Fe/H] at 
fixed [Z/H]. This makes the Fe I absorption lines in the 
continuum regions on either side of the index passband 
too strong, lowering the continuum measurement and 
consequently lowering the measured value of the Balmer 
absorption index in the models. This makes the observed 
Balmer index look stronger in comparison to the mod- 
els, and results in a younger age measurement. Using 
the correct [a/Fe] model fixes this problem, which has 
led Thomas, Maraston, & Korn (2004) to claim that a- 
enhancement effects cause younger ages to be measured 
in higher order Balmer absorption lines. However, the 
effect is dominated by Fe I lines. For models computed 
at fixed [Fe/H], the effect of the iron abundance on Hjp 
and HSp is accounted for, and does not depend on the 
a-enhancement of the model. Thus the claimed [a/Fe]- 
dependence of ages measured from Hy^ and HSf is in 
fact a dependence on [Fe/H] that is translated into an 
a-enhancement effect in models at constant [Z/H]. The 
S06 models, which are computed at fixed [Fe/H] rather 
than fixed [Z/H] do not have this entangled dependency. 

The S06 models produce consistent age measurements 
from H/3, Hjf, and KSp when they are used to analyze 
known SSPs such as the Galactic globular clusters 47 
Tuc, M5, and NGC 6528, and the open cluster M67 (Schi- 
avon 2007) which span a wide range of SSP parameters. 
Both 47 Tuc and M5 are metal-poor systems for which 
the a-enhancement effect on Hjf and HSf should be 
weak (Thomas, Maraston, & Korn 2004), and M67 is not 
significantly a-enhanced, thus these clusters are less rele- 
vant comparisons for high metallicity a-enhanced galax- 
ies. However, NGC 6528 is relatively old, has roughly 
solar [Fc/H] and is mildly a-enhanced, with age = 11 ±2 
Gyr, [Fe/H]= 0.0 ± 0.15, [Mg/Fe]= +0.10 ± 0.04 and 
[0/Fe]= +0.1 ± 0.05 (Zoccali et al. 2004, Barbuy et 
al. 2004, Carretta et al. 2001, summarized in Table 10 
of Schiavon 2007) and is thus an excellent SSP ana- 
log for early type galaxy data. The S06 best-fitting 
model for NGC 6528 has [Fe/H]= -0.2, [Mg/Fe]= +0.1, 
[0/Fc]= +0.15, and produces H/3, Hjf, and H5f ages 
of 13, 10, and 12 Gyr, respectively, which are relatively 
consistent within the uncertainties in the models. 

In the red sequence galaxies analyzed in Schiavon 
(2007), however, the models produce age measurements 
that decrease as bluer Balmer lines are used. Schiavon 
(2007) has shown that the most likely explanation for 
this phenomenon is that galaxies are not SSPs but con- 
tain stars with a range of ages. This effect has also been 
reported by Sanchez-Blazquez et al. (2006); they also 
interpret the differing age measurements from different 



Balmer lines as due to extended star formation histories. 

In Schiavon (2007), four possible explanations for 
this phenomenon are discussed: the presence of a sub- 
population of young to intermediate-age stars, a sig- 
nificant number of blue stragglers, the presence of a 
metal-poor old stellar population with a blue horizontal 
branch, or the effects of strongly non-solar abundance 
ratios which are not constrained by the available Lick 
indices. Among these scenarios, Schiavon (2007) con- 
cludes that a young-to-intermediate subpopulation is the 
most likely; the blue horizontal branch scenario cannot 
reproduce all three Balmer lines simultaneously, the frac- 
tion of blue stragglers required to reproduce the observed 
line strengths is 10-100 times larger than blue straggler 
populations observed in Galactic globular clusters, and 
the abundance ratios that reproduce the observed line 
strengths seem extreme: [O/Fe] = +1.0, or [Ti/Fe] = 
-0.7. 

In contrast, a small population of young-to- 
intermediate stars is sufficient to produce differences in 
Balmer line strengths similar to those observed. Schi- 
avon (2007) shows that a 11.2 Gyr SSP combined with 
a superimposed young population with age = 0.8 Gyr 
and comprising 0.5-1.0% of the total stellar mass pro- 
duces measured H/3 age ~ 5 Gyr and R5f age ~ 3.5 Gyr. 
These ages correspond to the measured H/3 and H<5^ ages 
of the cr = 70-120 km s" 1 strong LINER-likc emission 
stacked spectrum shown in Figure 12, with Ha/H/3 = 
3.27, which can therefore be modeled as such a two-burst 
system. This simple example serves to illustrate that 
measurement and modeling of multiple Balmer lines in 
a single system has the potential to reveal extended star 
formation histories. More sophisticated analysis of these 
galaxy spectra using composite stellar populations is be- 
yond the scope of this work but is clearly worthy of future 
investigation. 

It is noteworthy that almost all of the stacked spec- 
tra in our sample — quiescent and LINER-like, low a and 
high a — display this discrepancy in ages as measured by 
different Balmer lines. The only exception is the low- 
est a quiescent galaxy bin, for which H/3, Hy^, and HSp 
ages are relatively consistent at 8.4, 5.9, and 7.0 Gyr, 
respectively. 

One could postulate that the younger ages in bluer 
Balmer lines, instead of indicating a composite stellar 
population, were in fact caused by an inadequate treat- 
ment of a-enhancement in the S06 models, although we 
have argued above that Hyp and HSf are strongly af- 
fected by a-enhancement only in stellar population mod- 
els computed at fixed [Z/H]. A systematic variation of 
Hjf and H5f with [a/Fe] could produce the trend ob- 
served in the quiescent galaxies: that the higher a galax- 
ies (which have higher [Mg/Fe], as shown below in §5.3) 
show younger ages in H"fF and HSf than in H/3, but that 
all ages are nearly consistent in the lowest a galaxy bins 
which have nearly solar abundance ratios. However, it 
would be difficult to then explain the fact that the weak 
and strong LINER-like galaxies do show younger H'Jf 
and HSf ages at all a. In the next section, we will show 
that the weak and strong LINER-like galaxies have virtu- 
ally identical abundance ratios to the quiescent galaxies 
at the same a, thus it is difficult to invoke abundance 
ratio effects to fully explain the trends in Figure 12. 

If wc were to take the age trends in Figure 12 at face 
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value, the implication is that for quiescent galaxies, those 
with large a have a spread of stellar ages, while those 
a lower a are relatively close to SSPs. This would be 
in disagreement with compelling work by Thomas et al. 
(2005) who find more extended star formation histories 
in lower a early type galaxies than in higher a galaxies. It 
would also disagree with recent evidence from star form- 
ing galaxies at redshift z ~ 1 that lower mass galaxies 
have more extended star formation episodes (Noeske et 
al. 2007). However, younger ages in H7F and HSp are 
seen in Figure 12 for the galaxy bins with LINER- like 
emission, at all a. If the galaxies which have experi- 
enced relatively recent star formation episodes are pre- 
dominantly found in the weak and strong LINER-likc 
classes, the classification by emission line strength may 
have isolated the low a galaxies whose star formation 
ceased long ago and who may therefore be more likely to 
resemble SSPs. It is not immediately clear, however, why 
this would not also be the case for more massive galaxies. 
A more detailed investigation of complex model star for- 
mation histories would be necessary to further interpret 
the trends seen in Figure 12, which we defer to future 
work. 

Detailed analysis of the age differences between Balmer 
lines may in fact be pushing the data farther than is rea- 
sonable, given the uncertainties in the stellar population 
modeling process and the unknown contribution of oxy- 
gen. In addition, any conclusion that is strongly depen- 
dent on the a = 70-120 km s" 1 and a = 120-145 km s" 1 
bins should be regarded with caution, as those bins are 
known to be substantially incomplete (see §2.4). 

An important caveat that pertains to both this work 
and the Schiavon (2007) study is that each of these 
stacked spectra are composites of the spectra of several 
hundred galaxies. The differences in ages measured from 
various Balmer lines may therefore merely indicate that 
the galaxies included in each bin have a range of stellar 
population ages, rather than that each individual galaxy 
contains a composite stellar population. 

5.3. Abundances as a Function of a 

Iron abundances and elemental abundance ratios mea- 
sured with EZ_Ages are shown in Figure 14 as a func- 
tion of a. Plotted abundances are computed using a so- 
lar abundance isochrone with Ha/H/3 = 3.27; the abun- 
dances computed with Ha/H/3 = 4.1 or the a-enhanced 
isochrone (not shown) are very similar. As before, black, 
green, and red data points show results for quiescent, 
weak, and strong LINER-like galaxies respectively. The 
dotted line in each panel shows the solar value for refer- 
ence. Error bars are la errors computed in EZ_Ages by 
summing in quadrature the errors in each abundance due 
to individual Lick index measurement errors. The effect 
of errors in the fiducial ages and [Fe/H] are propagated 
through the abundance ratio estimates. Graves & Schi- 
avon (in preparation) will discuss the error calculations 
in greater detail. 

Unlike the mean, lightweighted SSP ages that were 
shown in Figure 12, the [Fe/H] and abundance ratios 
of the stacked galaxy spectra do not appear to vary with 
emission line strength; at fixed a, quiescent and LINER- 
like galaxies appear to have the same average abundance 
pattern within the errors. This means that the metal ab- 
sorption line variations with emission line strength that 



are apparent in Figure 6 are caused by variations in mean 
galaxy age, rather than by abundance differences. 

The chemical abundances do however depend on a. 
The solid black lines in Figure 14a-e are linear least 
squares fits of the abundances onto a. A single line is 
used to fit quiescent, weak, and strong LINER-like galax- 
ies simultaneously. For ease of comparison, the scale of 
the y-axis is the same in panels a-e. All abundances and 
abundance ratios increase with increasing a such that 
high a galaxies are both more metal-rich and more en- 
hanced in light elements with respect to Fe than are the 
low a galaxies. 

[Fe/H] increases with increasing a as [Fe/H] oc ct 45 , 
and is observed to range from about 1/2 solar at a ~ 100 
km s _1 to slightly sub-solar at a ~ 250 km s _1 . Be- 
cause early-type galaxies are known to have substantial 
mctallicity gradients (Carollo, Danziger, & Buson 1993, 
Ogando et al. 2006), the zero-point of the [Fe/H] relation 
is dependent upon the aperture used to extract the spec- 
tral data. The SDSS 1.5" spectral fibers typically cover 
0.35-1.1 effective radii and thus sample a large fraction 
of the galaxy light. Nuclear spectra of early- type galax- 
ies (e.g. those of Trager et al. 2000b and Thomas et al. 
2005) should appear substantially more metal-rich due 
to the metallicity gradients of the galaxies. 

The galaxies in this sample are Mg-cnhanced at all 
a, with the highest a galaxies being the most enhanced 
([Mg/Fe] oc ct 036 ). Both [C/Fe] and [N/Fe] increase 
strongly with c, from solar values in the low a galaxies 
to substantial C and N enhancement at high a. Finally, 
Ca appears to follow Fe in all galaxies. Table 7 gives the 
slope of the solid line fit for each elemental abundance 
or abundance ratio, along with the lcr uncertainty in the 
slope. 

The tendency of Ca to remain relatively constant with 
a, rather than to follow Mg and other a-elements as ex- 
pected if Ca production is dominated by Type II super- 
novae (Woosley & Weaver 1995), has been known for 
some time. It has been seen both in the Ca4227 index, 
which appears to be constant over a range of a (e.g., 
Vazdekis et al. 1997, Thomas, Maraston, & Bender 2003) 
and in the near-IR Ca 11 triplet, which may even anti- 
correlate with a (e.g., Saglia et al. 2002, Cenarro et al. 
2003, and references therein). Recently, Prochaska, Rose, 
& Schiavon (2005) have demonstrated that the Ca4227 
index is stronly affected by CN absorption in its blue 
continuum region and have suggested that redefining the 
Ca4227 index to avoid this region results in a Ca4227 
index that is less sensitive to other element abundances. 
Unfortunately, the Korn, Maraston, & Thomas (2005) 
index sensitivities and the S06 models do not include the 
redefined Ca4227r index. In principal, the abundance fit- 
ting process described in §4.2, which fits for [C/Fe] and 
[N/Fe] before using Ca4227 to fit for [Ca/Fe], should ac- 
count for the CN effect on Ca4227. Indeed, Schiavon 
(2007) uses the same abundance fitting method as in this 
paper and finds that [Ca/Fe] does increase with galaxy 
luminosity in a set of stacked SDSS spectra. However, 
the discrepancy between that result and our Figure 14 
appears to be an artifact due to imperfect emission in- 
fill correction in the Hf3 feature in the Schiavon (2007) 
analysis. Details of this are provided in appendix B. 

An interesting result from this Figure is that [C/Fe] 
and [N/Fe] show the strongest variation with a, with 
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considerably steeper slopes than [Mg/Fe]. The value of 
[Mg/Fe] is measured here using the Mg b Lick index, and 
has been used in past works (e.g., Trager et al. 2000a) 
to trace the enhancement of a-elements ([a/Fe]). If Mg 
is indeed a good proxy for the bulk of the a-elements, 
the relative strength of the abundance-cr relations shows 
that [C/Fe] and [N/Fe] vary more strongly with a than 
does [a/Fe]. 

It has been suggested in previous work that [C/Fe] 
and [N/Fe] are enhanced above solar ratios. Trager et al. 
(2000a) prefer a model for abundances in elliptical galax- 
ies that has C and N enhanced along with the a-elemcnts 
O, Ne, Na, Mg, Si, and S, but state that a similar model 
with solar C is difficult to distinguish from their pre- 
ferred model based on the data. Sanchez-Blazquez et 
al. (2006) prefer an abundance model in which the a- 
elements Na, Si, Ca, O, and Ti arc enhanced at higher 
a, with C enhanced but as a less steep function of a than 
the a-elements, with N enhanced more steeply than the 
a-elements, and with Mg enhanced most strongly of all. 
The results presented here differ from both Trager et al. 
(2000a) and Sanchez-Blazquez et al. (2006) in that we 
show [C/Fe] and [N/Fe] enhancement to increase more 
strongly with a than does [Mg/Fe]. All of these results 
are at odds with the C enhancement pattern found by 
Clemens et al. (2006), who find that [C/H] decreases with 
increasing a . It is clear that there is no consensus as to 
the enhancement (or depletion) pattern of these light el- 
ements. 

One difference in the modeling process is that the mod- 
els used in our analysis were computed at [O/Fe] = 0.0, 
whereas the other three works mentioned above have oxy- 
gen enhanced along with the a-elements. The dashed line 
in Figure 14a-e shows the effect of fixing [O/Fe] at +0.3 
during the abundance fitting process instead of [O/Fe] 
= 0.0. Both [Fe/H] and [Mg/Fe] appear unchanged and 
[Ca/Fe] shows only a small difference, but [C/Fe] and 
[N/Fe] increase at all a for the oxygen-enhanced mod- 
els. This is because oxygen has a negative effect on the 
C 2 4668 and CN index strengths used to compute [C/Fe] 
and [N/Fe]. In the presence of increased oxygen, most 
C will form CO rather than C2 or CN, requiring larger 
C and N abundances in order to produce the observed 
C2 and CN line strengths. Note that although there is 
a vertical offset, the slopes of the enhancement- a rela- 
tions do not change between the [O/Fe] = 0.0 and the 
[O/Fc] = +0.3 models. If [O/Fe] increases with a follow- 
ing [Mg/Fe] as expected if both Mg and O production 
are dominated by Type II supernovae (SNe), it is clear 
that the high a galaxies will be more enhanced above the 
[O/Fe] = 0.0 level than the low a galaxies, causing the 
slopes of the [C/Fe]-<r and [N/Fe]-<r relations to be even 
steeper. 

Many stellar population models are computed at fixed 
Z/H], whereas the S06 models are computed at fixed 
Fe/H]. To facilitate comparisons between models, Fig- 
ure 14f shows [Z/H] as a function of cr, for five different 
assumptions about [O/Fe] values. [Z/H] was computed 
using the solar elemental abundances of Grevesse, Noels, 
& Sauval (1996). In calculating [Z/H], Cr, Mn, Fe, Co, 
Ni, Cu, and Zn are treated as Fe-peak elements and as- 
sumed to scale with the measured [Fe/H], while Ne, Na, 
Mg, Si, S, and Ti are treated as a-elements and assumed 
to scale with the measured [Mg/Fe]. C, N, and Ca are 



allowed to vary separately with their measured abun- 
dance ratios. All other elements are assumed to have 
solar abundances, excepting O. The black lines in Figure 
14f show [Z/H] computed assuming constant O enhance- 
ment; [O/Fe] = 0.0, [O/Fc] = +0.3, and [O/Fc] = +0.5 
are shown as the solid, dashed, and dash-dot black lines, 
respectively. The blue lines show [Z/H] computed with 
[O/Fe] that varies with a. The solid blue line is com- 
puted assuming that O scales with the other a-elements 
([O/Fe] = [Mg/Fe]), while the dashed blue line has O 
scaling with Mg, but enhanced by a factor of two above 
Mg ([O/Fe] = [Mg/Fe] + 0.3). ' The slopes of the re- 
lations are given in Table 7. Unsurprisingly, a steeper 
relation results when O scales with Mg because [O/Fc] 
increases with a. Figure 14f highlights a problem with 
using models computed at fixed [Z/H]: total metallicity 
is dominated by O, which is highly uncertain and not 
measured. 

As discussed in Appendix A, Figure 15 shows that, at 
fixed cr, galaxies with fainter magnitudes have weaker 
metal line strengths than brighter galaxies but almost 
identical Balmer line strengths, suggesting that fainter 
galaxies are more metal-poor than their brighter coun- 
terparts at the same a. We have used the S06 models to 
convert the line strengths from Figure 15 into age and 
abundance estimates, following the same procedure de- 
scribed above for quiescent versus LINER- like galaxies. 
This modeling (not shown) confirms that fainter galaxies 
have lower values of [Fe/H] and older ages compared to 
brighter galaxies at the same cr, while their abundance 
ratios are nearly identical. 

5.4. Comparison with Previous Results 

To put the results presented here into the context of 
existing similar work, we have compiled results from four 
other SSP analyses of early type galaxies. Table 8 gives 
the slopes of the log age-log a, [Fe/H]-logcr, [Z/Hj-loga, 
and [a/Fe]-log(7 relations reported here and in previous 
work. Interestingly, although models claiming to accu- 
rately reproduce the effects of a-enhancement are rel- 
atively recent, the reported [a/Fe] values agree better 
than the ages or metallicities. There appears to be a 
consensus value of d[a/Fe]/d log a = 0.32 ± 0.04. 

All of the studies reviewed here find that typical galaxy 
ages increase with cr, although there is substantial vari- 
ation in the reported values of d log age/d log cr, ranging 
from 0.24 (Thomas et al. 2005) to 1.15 (Bernardi et al. 
2003b). Our value of d log age/d log a = 0.35 is toward 
the shallower end of the reported trends. It is noteworthy 
that the age trends shown in Figure 12 for LINER- like 
galaxies are steeper than those for quiescent galaxies. 

As for d[Z/H]/dlogcr, our value of 0.79 ± 0.05 com- 
puted assuming [0/Fc] = 0.0 (see Table 8) is in excel- 
lent agreement with the Trager et al. (2000b) value of 
0.76+0.13, and is a little steeper than the trends reported 
by Thomas et al. (2005), Nelan et al. (2005), Smith, 
Lucey, & Hudson (2007), and Bernardi et al. (2003b). 
Because this work analyzes stacked spectra, we can only 
quantify mean trends and thus the age-metallicity hy- 
pcrplanc of Trager et al. (2000b) is not visible in this 
work. 

6. DISCUSSION 
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We have shown that, on average, red sequence galax- 
ies with LINER-like emission have younger SSP ages 
than quiescent galaxies at the same a; strong LINER- 
like galaxies are 2-3.5 Gyr younger than their quiescent 
counterparts. This suggests that there is a connection 
between the star formation history of a galaxy and the 
presence or absence of LINER-like emission within it. 

Several caveats apply to this statement. Most impor- 
tantly, this study is based on stacked spectra, and there- 
fore analyzes only the mean properties of galaxies in each 
bin of emission properties and a. There may be a signif- 
icant spread in the ages of galaxies within each bin; not 
all LINER-like galaxies are necessarily younger than all 
quiescent galaxies at the same a. Thomas et al. (2005) 
show a spread in the ages of early-type galaxies at fixed 
a which they can match using Monte-Carlo simulations 
incorporating measurements errors and an intrinsic age 
spread of 20-25%. Trager et al. (2000b) also show a 
large spread in ages at fixed a. We have checked their 
reported ages as a function of [Om] emission line strength 
and confirm that the galaxies in their sample with > 2a 
[Om] emission line detections are generally younger than 
their quiescent galaxies at similar a. It is possible that 
only a subset of LINER-like galaxies have younger ages, 
but these must then make up a significant fraction of the 
LINER-like galaxies in order to substantially lower the 
mean LINER-like galaxy ages to match observations. 

It is also clear from Figure 12 that there cannot be 
a one-to-one correlation between SSP age and LINER- 
like emission strength because the mean H/3 age of the 
highcst-CT LINER-like galaxies is older than the mean 
H/3 age of lower-cr quiescent galaxies; specifically, if 
LINER-like emission existed in all galaxies with ages be- 
low some threshold value, all LINER-like galaxies would 
be younger than all quiescent galaxies, regardless of a, 
which is not seen. However, this trend disappears us- 
ing H6p, where all mean ages of LINER-like galaxies are 
younger than all the HSf mean ages of quiescent galax- 
ies, even though this is not true of H/3. If, as suggested 
in Schiavon (2007), higher order Balmer lines are more 
sensitive to the presence of very young stars, this may 
imply that LINER-like galaxies have a larger fraction of 
relatively younger stars than quiescent galaxies, that is 
the LINER-like galaxies have had more recent star for- 
mation than their quiescent counterparts. 

This does not mean that the LINER-like galaxies are 
post-starburst systems, galaxies which have undergone a 
recent epoch of powerful star formation but are no longer 
forming stars. Spectra of post-starburst galaxies show 
strong higher-order Balmer absorption lines characteris- 
tic of A stars. They are often modeled by the linear com- 
bination of an A star spectrum and a K star spectrum. 
The ratio of the linear coefficients of this two-component 
model (A/K) is used to make the post-starburst classifi- 
cation. Only 1.1% of the LINER-like galaxies have A/K 
ratios from Yan et al. (2006) consistent with identifica- 
tion as post-starburst galaxies (A/K > 0.25), suggesting 
that the LINER-like galaxy population in general is more 
than ~ 0.5-1 Gyr beyond any major recent episode of 
star formation (Quintero et al. 2004). However, LINER- 
like galaxies might be "post post-starburst" systems seen 
> 1 Gyr after the post-starburst phase, possibly a con- 
tinuation of the post-starburst AGN population seen in 
Kauffmann et al. (2003) to older ages and lower ioniza- 



tion AGN. It is also possible that the LINER-like galax- 
ies are undergoing low-level star formation which is con- 
tributing to the optical emission lines. 

The interpretation of the age difference between quies- 
cent and LINER-like galaxies is not straightforward, due 
to the lack of consensus as to the ionization mechanism 
driving LINER-like emission; it is in fact likely that LIN- 
ERs and LINER-like galaxies arc a heterogeneous pop- 
ulation (Filippenko 2003). Possible mechanisms for pro- 
ducing LINER-like emission are: 

1. Photoionization by an AGN continuum (e.g., Fer- 
land & Netzer 1983). 

2. Photoionization by newly-formed O stars (e.g., 
Shields 1992). 

3. Shock ionization (e.g., Dopita & Sutherland 1995). 
The shocks could be driven by starburst winds, gas 
accretion from a satellite galaxy or the surrounding 
inter-galactic medium, episodic AGN feedback, or 
by galaxy interactions and/or mergers. 

4. Cooling flows (e.g., Heckman 1981). 

5. Photoionization at early times (< 500 Myr, e.g., 
Taniguchi, Shioya, & Murayama 2000) or at late 
times (1-13 Gyr, e.g., Binctte et al. 1994) by post- 
AGB stars. 

Of these possibilities, photoionization or shock ioniza- 
tion by young stars and photoionization by post-AGB 
stars are the most directly associated with the stellar 
population of the galaxy, and would thus give a poten- 
tially straight-forward connection between the presence 
of LINER-like emission and the star formation history 
of a galaxy. In particular, if the observed LINER-like 
emission in many of these galaxies is powered by young 
stars, it would not be at all surprising to discover that 
these galaxies have younger mean light-weighted ages. 

However, there are several indications that this expla- 
nation is too simple. First, Yan et al. (2006) show that 
the LINER-like red sequence galaxies in their analysis 
have higher observed [On] /H/3 ratios than either star- 
forming galaxies or AGN (i.e. galaxies with Seyfert- 
like emission line ratios), and thus LINER-like emission 
cannot be produced by a superposition of Hn regions 
and AGN activity. This does not preclude the LINER- 
like emission being caused by star formation through a 
mechanism other than Hn region emission, however Hn 
regions produce strong H/3 emission and their presence 
would tend to weaken the [On] /H/3 ratios of the galaxy. 

There are galaxies with LINER-like emission that do 
show indications of containing starbursts. These are the 
IR-luminous galaxies described in Vcilleux et al. (1995) 
and Sturm et al. (2006). These galaxies have mid-IR 
SEDs and IR fine structure emission line ratios similar to 
those observed in starburst galaxies (Sturm et al. 2006). 
In addition, they have highly reddened but intrinsically 
blue optical continuum colors and low Mg b EWs char- 
acteristic of relatively young stellar populations in dusty 
galaxies (Veilleux et al. 1995). We do not have mid- 
IR data for the SDSS galaxies to determine whether the 
galaxies present in our sample arc indeed IR-luminous, 
but the optical continuum colors, internal reddening (de- 
termined by observed Ha/H/3 emission line ratios), and 
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Mg b EWs of the IR-luminous LINER-like galaxies are 
vastly different from those measured in the galaxy sam- 
ple presented here. We therefore think it highly unlikely 
that these IR-luminous LINER-like galaxies contribute 
significantly to the LINER-like red sequence galaxies in 
this sample. 

In the scenario where the bulk of the LINER-like emis- 
sion on the red sequence is produced by post-AGB star 
photoionization, a correlation between LINER-like emis- 
sion and stellar population age would imply that either 
the amount of ionizing photons produced by the old stel- 
lar population or the supply of gas (from stellar mass 
loss) would depend upon the age of the stars in the 
galaxy, which is predicted by some models (e.g., Mathews 
1989). 

The other possible LINER mechanisms mentioned 
above (photoionization by an AGN, shock ionization due 
to gas accretion, galaxy interactions, or merger activity, 
and cooling flows) are not directly related to the stellar 
populations of the galaxies. To explain the correlation 
between SSP age and emission properties with any of 
these mechanisms requires a connection between the stel- 
lar population age of the galaxy and some quantity exter- 
nal to the galaxy (accretion, interaction, cooling flows) or 
internal to the galaxy but separate from the stellar pop- 
ulation (AGN activity). This would also be true of a sce- 
nario in which the post-AGB stars in the galaxy provide 
the ionizing photons but where the emission strength is 
determined by the supply of gas to be ionized, with the 
gas coming from an external source (rather than from 
stellar mass loss). 

These "external" possibilities have interesting implica- 
tions for galaxy evolution, in particular the possibility 
of ionization by low luminosity AGN. This mechanism is 
favored by Filippenko (2003) for centrally-concentrated 
LINER emission, and has a tantalizing resonance with 
recent work suggesting that AGN activity is involved in 
quenching star formation in galaxies (e.g., Kauffmann et 
al. 2003, Hopkins et al. 2005, Croton et al. 2006, and oth- 
ers). However, there is some evidence that the LINER- 
like emission in many if not most early-type galaxies is 
extended throughout the galaxy, rather than being cen- 
trally concentrated (e.g., Goudfrooij et al. 1994, Sarzi 
et al. 2006), suggesting that photoionization by an AGN 
may not be the dominant mechanism driving LINER-like 
emission in most early-type galaxies. 

Unfortunately, it is outside the scope of this work to 
speculate as to which mechanism (or mechanisms) dom- 
inate the observed LINER-like emission in red sequence 
galaxies, making it difficult to assess the physical cause of 
our observed result. It would be particularly interesting 
to perform an analysis similar to that presented here on 
a sample of galaxies with spatially resolved spectra, so 
that true nuclear LINERs could be analyzed separately 
from galaxies with extended LINER-like emission. 

7. CONCLUSIONS 

This paper has presented an analysis of stacked spec- 
tra for a well-defined set of SDSS red sequence galaxies. 
Galaxies observed to harbor LINER-like emission are an- 
alyzed in parallel with quiescent galaxies, while galax- 
ies with emission line ratios characteristic of Hn regions, 
Seyferts, or Transition Objects have been excluded. Lick 
index absorption line strengths have been measured in 



all the stacked spectra and have been used to study the 
stellar population properties of red sequence galaxies as 
a function of galaxy velocity dispersion and emission line 
properties. The main results of this analysis can be sum- 
marized as follows: 

1. Red sequence galaxies with LINER-like emission 
(roughly one-third of the total red sequence) are 
systematically younger than their quiescent coun- 
terparts at similar a by several Gyr. Despite this 
age difference they appear to have metallicities and 
abundance patterns similar to quiescent galaxies. 
This result is robust to the details of the emis- 
sion infill correction, the enhancement pattern of 
the isochrone used to construct comparison mod- 
els, and the (unknown) oxygen abundance of the 
galaxies. 

2. Many properties of red sequence galaxies' stellar 
populations vary with galaxy velocity dispersion, 
including SSP age, [Fe/H], and abundance ratios. 
Galaxies with higher a are typically older and more 
metal-rich than lower a galaxies, and are increas- 
ingly enhanced in Mg, C, and N. Ca appears to 
scale with Fe at all a. 

3. [C/Fe] and [N/Fe] vary strongly with a; the slopes 
of the [C/Fe]-log a and [N/Fc]-log a relations are 
steeper than the [Mg/Fe]-log a relation. 

4. The SSP ages measured in stacked spectra from 
the bluer Balmer lines (H7^ and H<5p) are younger 
than those measured using H/3. This effect is not 
present when the S06 SSP models are used to ana- 
lyze known simple stellar populations. If the inter- 
pretation presented in Schiavon (2007) is correct, 
the discrepant age measurements observed in this 
analysis indicate the presence of composite stellar 
populations. From the work presented here, it is 
not possible to tell whether the individual galaxies 
that go into each stacked spectrum are themselves 
composite stellar populations with similar star for- 
mation histories, or whether the individual galax- 
ies are true SSPs of different ages which produce 
a composite population when they are combined 
together. 

The correlation between the presence of LINER-like 
emission and intrinsically younger stellar populations 
demonstrates that there is a link between the mecha- 
nism driving the emission and the star formation history 
of the galaxy. If the LINER-like emission in red sequence 
galaxies is primarily driven by AGN activity, this may be 
an indication that AGN do indeed play a role in ending 
star formation activity. Alternatively, if the LINER-like 
emission in most of these galaxies is controlled entirely 
by the properties of the stellar population itself (e.g. gas 
lost from evolved stellar envelopes which is ionized by 
UV flux from post-AGB stars), the correlation with stel- 
lar population age may provide important constraints on 
models for the UV flux from low mass stars or for the 
rate of mass loss in a stellar population over time. 
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APPENDIX 

A. CORRECTING FOR INCOMPLETENESS AT LOW a 

In §2.4, we showed that the magnitude limit of the SDSS makes the galaxy sample presented in this paper significantly 
incomplete in the lower a bins, with the lowest a bin missing as much as 51% of the galaxies it should contain. This 
causes the low a bins to be increasingly biased toward the brightest galaxies in that bin. 

To understand the effect of this bias on the measured line strengths, it is necessary to characterize the variation 
of mean absorption line strengths as a function of galaxy magnitude at fixed a. We therefore use the total sample 
of 10,284 galaxies with no Ha or [On] emission detected at the 2a level to construct high S/N stacked spectra as 
a function of a and absolute magnitude 0A M r . As in the selection of the sample of 2000 quiescent galaxies, we 
exclude objects with possible low-level emission and maintain a roughly symmetric distribution about EW([On]) = 
by requiring -3.12 A< EW([On]) < 3.12 A (see §2.3). This leaves us with a sample of roughly 9400 galaxies. 

These galaxies are divided into six bins in a, as before. Each of these six bins is then divided into three bins in 
absolute magnitude, with equal numbers of galaxies in each magnitude bin. The median magnitudes of galaxies in the 
bright bins are ~ 0.5 mag brighter than those in the middle bins at the same a, which are in turn ~ 0.5 mag brighter 
than those in the faint bins. Galaxies are coadded within each bin (using the process described in §2.3), and the Lick 
indices and D n 4000 are measured in each of the cr-magnitude binned spectra. 

Measured line strengths are shown as a function of a in Figure 15. In each a range, the brightest, middle, and 
faintest magnitude bins are shown in red, green, and blue, respectively. As in Figure 6, error bars include statistical 
errors only, and the lines show a least squares fit of index strength onto a separately for the brightest, middle, and 
faintest galaxies at each a. For ease of comparison, the axis scales in each index plot are the same as in Figure 6. 

From Figure 15, it is clear that the stellar populations of red sequence galaxies do vary with absolute magnitude 
at fixed a. The iron lines (Figure 15b) in particular are considerably stronger in bright galaxies than in faint ones. 
However, the Balmer lines and G4300 (Figure 15a) show little variation, while other lines (Figure 15cd) show some 
slight variation with magnitude. This suggests that at fixed a, brighter galaxies have higher metallicities and slightly 
younger ages compared to fainter galaxies. 

The low a bins of the quiescent and LINER-like galaxy sample are biased toward brighter objects and therefore, 
because of the variation in line strengths with magnitude, toward stronger absorption lines for all lines shown in Figure 
15b-d. A correction for this effect is computed as follows: we assume that the characteristics of the galaxies missing 
from each a bin are similar to the faint galaxies that are included in the bin, and that the mean measured properties 
of the bin are similar to those of the middle magnitude bin. For each measured index, the difference between the 
middle magnitude bin and the faint bin at fixed a is computed by subtracting the linear fit to the middle magnitude 
bins (green lines in Figure 15) from the faint bins (blue lines). The line strength differential at each a is multiplied by 
the corresponding incompleteness fraction (see Table 1), so that the a = 70-120 km s _1 bin has the largest correction 
because it is most incomplete, while the a = 220-300 bin has no correction because it is complete by assumption (see 
§2.4). This scaled correction is then subtracted from the measured line strengths in the quiescent, weak, and strong 
emission bins at each value of a. The scaled corrections for each index are given in Table 6 for the five lowest a bins 
(the highest a bin is assumed to be complete, so no correction is applied) . Implicit in this treatment is that the relative 
line strength differences with magnitude are the same for LINER-like galaxies as for quiescents. 

The corrected line strengths are shown in Figure 6 in the right-hand column of each panel. Comparing the corrected 
with the uncorrected (left-hand column) values, it is clear that the corrections are small and are negligible in the 
Balmer lines and G4300. In the case of the metal lines (Figure 6b-d), they operate in the sense of slightly enchanting 
trends in line strength with a. They have no effect on the relative line strengths of galaxy bins with differing emission 
properties because the same correction was applied to all bins at the same a. 

B. COMPARISON WITH RESULT FROM SCHIAVON 2007 

Although the results presented in this work are generally consistent with the analysis of stacked SDSS spectra in 
Schiavon (2007), there are a number of differences worthy of deeper investigation. Schiavon (2007) finds no trend in 
age with galaxy luminosity when using H/3 to compute mean ages, but finds that galaxy age increases with increasing 
luminosity when ages are computed using H<5p. In contrast, we find age increasing with a regardless of which Balmer 
line is used in the fitting process. In addition, Schiavon (2007) finds that [Ca/Fe] increases with galaxy luminosity, 
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whereas we find practically no trend in [Ca/Fe] with a. It is difficult to directly interpret these differences, as the 
Schiavon (2007) results are presented as a function of galaxy absolute magnitude, rather than a. 

In order to make a more direct comparison between the results of Schiavon (2007) and those presented here, we 
have taken the same galaxy sample used in our main analysis and divided the galaxies into bins based on absolute 
magnitude, rather than a. As before, galaxies with weak and strong LINER-like emission are stacked separately 
from those with no detectable emission lines. For this process, we have used the galaxy absolute r-band magnitude, 
K-corrected to z=0.0 in order to match the Schiavon (2007) data which made use of stacked spectra from Eisenstein et 
al. (2003). We then measure Lick indices and fit for stellar population mean ages, abundances, and abundance ratios 
as described in §3 and §4 above. 

Figure 16 shows the results of the stellar population modeling for our luminosity-binned sample (black stars), as 
compared to the results from Schiavon (2007) (grey triangles). Values of [Fe/H] appear slightly higher in the Schiavon 
(2007) results than in our luminosity-binned sample, by a few hundredths of a dex. The typical values for the abundance 
ratios are comparable between the two data sets, but the Schiavon (2007) values show consistently stronger trends 
with magnitude than do our data. Again, the discrepancies are all within a few hundredths of a dex. The largest 
discrepancies between the data sets are in the age measurements. Here our data show slightly older ages for more 
luminous galaxies, while the Schiavon (2007) results show slightly younger ages for more luminous galaxies. 

Because both the Eisenstein et al. (2003) stacked spectra and the stacked spectra in our sample were acquired as 
part of the same survey and have been analyzed using the same method, any differences in the results should be due 
to differences in the sample selection process. The sample selection criteria are in fact quite different between the 
two samples. Both samples are taken from the SDSS MAIN galaxy sample. The Eisenstein et al. (2003) spectra 
used in Schiavon (2007) are chosen using two morphological criteria: that the galaxies have concentrated light profiles 
(determined from the SDSS concentration parameter c, see Eisenstein et al. 2003 for details) and by requiring that a de 
Vaucouleurs fit to the light profile be 20 times more likely than an exponential fit. Most of these galaxies lie on the red 
sequence; an additional color criterion is applied to exclude color outliers, removing about 5% of the morphologically 
selected galaxies from the sample. 

It is important to highlight here that no spectral criteria are applied. Many of the galaxies in the Eisenstein et al. 
(2003) stacked spectra may therefore include emission lines from ionized gas. In fact, the principle component analysis 
in Eisenstein et al. (2003) shows that there is a significant component of emission in the galaxies in their sample. 
Yan et al. (2006) have shown that red sequence galaixes (chosen by a color criterion only, without any morphological 
selection) consist of about 48% quiescent galaxies without emission lines, 29% galaxies with LINER-like emission, 
and 23% galaxies with emission characteristic of Seyferts or Hn regions. The morphological selection applied by 
Eisenstein et al. (2003) may remove some fraction of the Seyferts and star-forming galaxies from the red sequence, but 
a contaminating population of galaxies with emission due to low-level star formation or Seyfert activity likely exists 
at the roughly 10% level. 

Before measured Lick indices and analyzing the stellar populations of the stacked Eisenstein et al. (2003) spectra, 
Schiavon (2007) corrects for emission infill in the stacked spectra by measuring the strength of the [On]A3727 line and 
applying the conversion from [On] EW to Ha EW determined in Yan et al. (2006). A standard Balmer decrement is 
then used to correct H/3, H"fF, and H<5p . However, this conversion is appropriate only for galaxies with LINER-like 
emission line ratios. The Eisenstein et al. (2003) spectra may contain some fraction of Seyferts and low-level star 
formation. For these galaxies, the strength of Ha emission is much higher for a given [On] EW, thus the emission infill 
correction performed by Schiavon (2007) may undercorrect for infill in the Balmer lines, to the extent that galaxies 
in the sample have non-LINER-likc emission line ratios. As discussed in §3.3, H/3 will be much more sensitive to the 
effects of infill than will the bluer Balmer lines. 

Figure 17 shows Lick index measurements as a function of galaxy absolute magnitude for our sample (black stars) 
and for the Eisenstein et al. (2003) stacked spectra (grey triangles). Only Lick indices that are used in the abundance 
fitting process are shown. The Balmer lines (left panel) differ significantly between the two data sets, while the metal 
lines Fe5270, Fe5335, Mg b, C24668, CN2, and Ca4227 are relatively consistent. In our luminosity-binned data, all 
three Balmer line strengths decrease as galaxy luminosity increases. The same trend is seen in HSp and H7^ for the 
Eisenstein et al. (2003) data, although the line strengths are offset toward stronger absorption compared to our data. 
This suggests that the Eisenstein et al. (2003) sample of galaxies includes a larger fraction of intermediate-age galaxies 
than does our own. 

This offset in mean stellar population age between the samples is not unexpected. The galaxies included in our 
luminosity-binned sample, as in the c-binned sample used in the main analysis presented above, have been selected 
with a fairly stringent color cut, which excludes the blue side of the red sequence (see Figure 1). No such stringent 
cut has been applied to the Eisenstein et al. (2003) galaxies. That sample therefore contains a sizable fractin of 
galaxies which are as much as 0.1 mag bluer in g — r than our galaxy sample and therefore may include a significant 
sub-population of intermediate-age galaxies which are excluded from our sample. 

Also, as discussed earlier, the Eisenstein et al. (2003) data likely include a modest fraction of galaxies with low-level 
star formation because galaxies have not been excluded on the basis of emission lines. Including these galaxies would 
bias the sample toward younger mean ages. It is true that the star-forming galaxies also have stronger Balmer emission 
infill, which we have argued is likely to be under-corrected. However, the effect on the blue Balmer lines is small both 
in terms of the total amount of infill and the amount of correction applied, thus the inaccuracies in accounting for the 
infill in RSf and H7^ should not be large. 

This is not true for H/3, where the emission infill is substantial and therefore the effect of the under-correction could 
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be significant. Indeed, H/3 behaves differently in the Eisenstein et al. (2003) stacked spectra than in our luminosity- 
binned spectra; H/3 remains roughly constant over a range of galaxy absolute magnitude in the Eisenstein et al. (2003) 
data whereas H/3 strengths decline for brighter galaxies in our sample. The contant values of H/3 are also qualitatively 
different from the HSf and Hjp indices in both the Eisenstein et al. (2003) data and in our data, where Balmer 
absorption declines with increasing galaxy luminosity everywhere. It seems likely that the H/3 values measured in the 
Eisenstein et al. (2003) data suffer from an improper infill correction which has dramatic effects on H/3 but less so on 
H5f and H-Jf, and none at all on the metal lines. 

The metal lines in fact appear very similar between the two different data sets. There is some indication that 
metal absorption strengths are slightly lower in the Eisenstein et al. (2003) data, which would be consistent with the 
supposition that those data include a larger fraction of intermediate-age galaxies than do our data. The slopes of the 
relations between line strengths and galaxy magnitude are in good agreement between the two samples. 

We therefore conclude that the slight differences in the ages and abundance patterns reported by Schiavon (2007) 
as compared to the results presented in this paper are due to differing trends in the H/3 line strengths that are used 
to compute the fiducial galaxy ages (see §4), rather than genuine differences in abundance patterns. These differing 
trends in measured H/3 line strengths are due to differences in sample selection coupled with less accurate emission infill 
corrections applied to the Eisenstein et al. (2003) data in Schiavon (2007). In order for this to explain the observed 
flat relationship between H/3 absorption strength and galaxy luminosity, the fainter galaxies would have to be more 
affected by under-corrected emission infill (in other words, the fainter galaxies would have to have a larger fraction 
of galaxies with emission having non-LINER-like line ratios). It seems reasonable that the fainter Eisenstein et al. 
(2003) galaxies, which are observed to have much stronger H<5p absorption than the brighter galaxies, are more likely 
to contain residual star formation and weak Hn region emission than the brighter galaxies. 

It should be noted that the differences in abundances shown in Figure 16 between the two data sets are never 
more than 0.05 dex. Although different sample selection and inaccurate emission infill corrections pose problems for 
precision stellar population determinations, there does not seem to be a catastrophic discrepancy created by moderate 
differences in the input data. 
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204.1 


-0.11 






weak 


374 


561.9 


203.9 


3.25 






strong 


347 


506.7 


203.2 


9.06 


220-300 




quiescent 


284 


526.4 


245.8 


-0.27 






weak 


289 


570.9 


247.1 


3.14 






strong 


326 


550.6 


243.6 


9.67 



TABLE 2 

Indices in Stacked Spectra: Balmer Indices* 



H/3 H 7F H<5 F 
Name (A) (A) (A)_ 



Infill Correction:** 


Ho97 


Yan06 


none 


Ho97 


Yan06 


none 


Ho97 


Yan06 


none 


70-120 


/ 


quiescent 


1.88 


1.88 


1 


.88 


-0.90 


-0.90 


-0.90 


0.58 


0.58 


0.58 


70-120 


/ 


weak 


1.95 


1.90 


1 


.69 


-0.45 


-0.48 


-0.61 


1.04 


1.02 


0.95 


70-120 


/ 


strong 


2.15 


2.03 


1 


.47 


-0.03 


-0.12 


-0.45 


1.39 


1.33 


1.14 


120-145 


/ 


quiescent 


1.79 


1.79 


1 


.79 


-1.06 


-1.06 


-1.06 


0.53 


0.53 


0.53 


120-145 


/ 


weak 


1.93 


1.89 


1 


.67 


-0.68 


-0.72 


-0.85 


0.93 


0.91 


0.83 


120-145 


/ 


strong 


2.04 


1.93 


1 


.43 


-0.47 


-0.55 


-0.85 


1.03 


0.98 


0.80 


145-165 


/ 


quiescent 


1.78 


1.78 


1 


.78 


-1.24 


-1.24 


-1.24 


0.49 


0.49 


0.49 


145-165 


/ 


weak 


1.87 


1.82 


1 


.62 


-0.89 


-0.92 


-1.04 


0.79 


0.77 


0.69 


145-165 


/ 


strong 


1.98 


1.85 


1 


.32 


-0.73 


-0.82 


-1.15 


0.85 


0.80 


0.61 


165-190 


/ 


quiescent 


1.73 


1.73 


1 


.73 


-1.33 


-1.33 


-1.33 


0.47 


0.47 


0.47 


165-190 


/ 


weak 


1.83 


1.78 


1 


.58 


-1.07 


-1.11 


-1.23 


0.70 


0.68 


0.60 


165-190 


/ 


strong 


1.93 


1.81 


1 


.29 


-0.98 


-1.07 


-1.39 


0.73 


0.68 


0.49 


190-220 


/ 


quiescent 


1.70 


1.70 


1 


.70 


-1.45 


-1.45 


-1.45 


0.32 


0.32 


0.32 


190-220 


/ 


weak 


1.73 


1.68 


1 


.49 


-1.29 


-1.32 


-1.44 


0.51 


0.49 


0.42 


190-220 


/ 


strong 


1.83 


1.70 


1 


.16 


-1.19 


-1.28 


-1.63 


0.63 


0.57 


0.37 


220-300 


/ 


quiescent 


1.63 


1.63 


1 


.63 


-1.61 


-1.61 


-1.61 


0.22 


0.22 


0.22 


220-300 


/ 


weak 


1.63 


1.59 


1 


.40 


-1.51 


-1.54 


-1.66 


0.35 


0.33 


0.26 


220-300 


/ 


strong 


1.72 


1.58 


1 


.01 


-1.41 


-1.51 


-1.88 


0.47 


0.40 


0.18 



Index values arc as measured in the stacked spectra. They do not include the corrections for incompleteness 
described in Appendix A Balmer index measurements are shown for three different choices of emission infill 
correction: (1) assuming Ha/H/3 — 3.27 as in Ho, Filippcnko, Sz Sargent (1997), (2) assuming Ha/H/3 — 4.1 
as in Yan et al. (2006), and (3) with no infill correction. See §2.2 for details. 



TABLE 3 

Indices in Stacked Spectra: D„4000, Fe, and Mg Indices 1 











D n 4000 


Fe4383 


Fe4531 


Fc5270 


Fe5335 


Fe5406 


Mgi 


Mg 2 


Mg b 




Name 




(A) 


(A) 


(A) 


(A) 


(A) 


(mag) 


(mag) 


(A) 


t u- 


1 on 
-IzU 


/ 


quiescent 


1 07 


q no 
o.yz 


o Qn 
z.oy 


o aA 
z.04 


O QO 


1 

l.OD 


n nQn 
U.UoU 


n one; 


Q A Q 


f U" 


-120 


/ 


weak 


1 Q1 

l.ol 


o on 

o.oo 


O Q7 
Z.O / 


o c;n 
z.oU 


o ot; 

z.zo 


1 ^O 
l.OZ 


n n7^ 
U.U / 


n 1 nc; 

u. iyo 


o.o5 


*7n 
( u- 


-120 


/ 

/ 


strong 


L. to 


q e;o 
o.oz 


O QQ 

z.oo 


o An 


O 1 Q 

z. lo 


1 A ft 
1.4D 


n nfio 

u.uoy 


n 1 qo 
U. loz 


Q 1 Q 

o.lo 


i on 
IzU- 


-145 


/ 


quiescent 


i.yi 


a 1 n 

4. iy 


q nfi 


o 7n 
z. i u 


O QO 

z.oy 


1 dO 


n nm 

u.uy 1 


n ooq 


o. / D 


i on 


-145 


/ 


weak 


1 or; 
l.oo 


q on 

o.yy 


o nn 

z.yy 


o 7n 
z. / U 


o /in 
Z.4U 


1.D4 


n no c; 


n oi q 
U.zlo 


Q fiO 


120- 


-145 


/ 


strong 


1.82 


3.92 


2.99 


2.65 


2.29 


1.50 


0.083 


0.206 


3.58 


145- 


-165 


/ 


quiescent 


1.94 


4.28 


3.03 


2.71 


2.46 


1.64 


0.100 


0.237 


3.91 


145- 


-165 


/ 


weak 


1.89 


4.16 


3.04 


2.70 


2.45 


1.63 


0.094 


0.225 


3.76 


145- 


-165 


/ 


strong 


1.87 


4.13 


3.03 


2.69 


2.42 


1.60 


0.093 


0.222 


3.82 


165- 


-190 


/ 


quiescent 


1.96 


4.25 


3.13 


2.74 


2.53 


1.67 


0.105 


0.244 


4.10 


165- 


-190 


/ 


weak 


1.90 


4.25 


3.10 


2.75 


2.47 


1.65 


0.100 


0.237 


4.02 


165- 


-190 


/ 


strong 


1.92 


4.31 


3.07 


2.76 


2.46 


1.66 


0.103 


0.239 


4.05 


190- 


-220 


/ 


quiescent 


1.99 


4.35 


3.16 


2.81 


2.52 


1.73 


0.114 


0.254 


4.16 


190- 


-220 


/ 


weak 


1.99 


4.42 


3.19 


2.79 


2.57 


1.70 


0.112 


0.256 


4.27 


190- 


-220 


/ 


strong 


1.98 


4.40 


3.15 


2.70 


2.50 


1.67 


0.112 


0.252 


4.27 


220- 


-300 


/ 


quiescent 


2.06 


4.48 


3.28 


2.81 


2.65 


1.70 


0.126 


0.276 


4.55 


220- 


-300 


/ 


weak 


2.01 


4.54 


3.23 


2.85 


2.64 


1.77 


0.125 


0.270 


4.45 


220- 


-300 


/ 


strong 


1.99 


4.45 


3.23 


2.84 


2.58 


1.75 


0.126 


0.271 


4.56 



Index values are as measured in the stacked spectra. They do not include the corrections for incompleteness 
described in Appendix A 



TABLE 4 

Indices in Stacked Spectra: CN, CH, C2, Ca, Na, and TiO Indices* 









CNi 


CN 2 


G4300 


C 2 4668 


Ca4227 


Ca4455 


NaD 


TiOi 


Ti0 2 


Name 


(mag) 


(mag) 


(A) 


(A) 


(A) 


(A) 


(A) 


(mag) 


(mag) 


70-120 


/ 


quiescent 


0.0263 


0.0507 


4.83 


5.40 


0.95 


1.12 


2.61 


0.0248 


0.0612 


70-120 


/ 


weak 


0.0139 


0.0398 


4.56 


5.06 


0.83 


1.02 


2.64 


0.0246 


0.0604 


70-120 


/ 


strong 


-0.0051 


0.0207 


4.19 


4.63 


0.86 


0.97 


2.77 


0.0246 


0.0573 


120-145 


/ 


quiescent 


0.0386 


0.0659 


4.99 


5.70 


0.99 


1.16 


3.08 


0.0271 


0.0643 


120-145 


/ 


weak 


0.0261 


0.0528 


4.74 


5.62 


0.93 


1.07 


2.98 


0.0234 


0.0627 


120-145 


/ 


strong 


0.0184 


0.0451 


4.59 


5.26 


0.91 


1.07 


3.24 


0.0238 


0.0614 


145-165 


/ 


quiescent 


0.0468 


0.0742 


5.05 


6.02 


1.00 


1.21 


3.18 


0.0280 


0.0677 


145-165 


/ 


weak 


0.0377 


0.0639 


4.82 


5.91 


0.94 


1.16 


3.41 


0.0241 


0.0655 


145-165 


/ 


strong 


0.0357 


0.0626 


4.85 


5.68 


0.93 


1.10 


3.40 


0.0249 


0.0650 


165-190 


/ 


quiescent 


0.0559 


0.0859 


5.17 


6.13 


0.97 


1.21 


3.18 


0.0254 


0.0689 


165-190 


/ 


weak 


0.0487 


0.0769 


4.96 


6.22 


0.96 


1.18 


3.48 


0.0256 


0.0675 


165-190 


/ 


strong 


0.0478 


0.0763 


4.97 


6.06 


0.95 


1.14 


3.50 


0.0265 


0.0689 


190-220 


/ 


quiescent 


0.0694 


0.1025 


5.22 


6.55 


1.02 


1.23 


3.62 


0.0286 


0.0713 


190-220 


/ 


weak 


0.0640 


0.0938 


5.11 


6.63 


1.02 


1.21 


3.83 


0.0280 


0.0711 


190-220 


/ 


strong 


0.0629 


0.0929 


5.14 


6.43 


1.01 


1.21 


3.96 


0.0272 


0.0712 


220-300 


/ 


quiescent 


0.0896 


0.1233 


5.32 


6.91 


1.08 


1.32 


4.16 


0.0321 


0.0762 


220-300 


/ 


weak 


0.0842 


0.1168 


5.20 


6.99 


1.03 


1.28 


4.47 


0.0305 


0.0754 


220-300 


/ 


strong 


0.0823 


0.1134 


5.20 


6.87 


1.03 


1.26 


4.41 


0.0290 


0.0744 



Index values are as measured in the stacked spectra. They do not include the corrections for incompleteness 
described in Appendix A 



TABLE 5 

Direction of Line Strength Change in SSP Models 



Balmcr lines D„4000, G4300 Metal lines 



Agef ^ ii ft 

[Fc/H] ft 41- ft It 



24 



Graves et al. 



TABLE 6 

Incompleteness Corrections 



a bin a bin a bin <r bin <j bin 



Index 


70-120 km s' 1 


120-145 km s" 1 


145-165 km s" 1 


165-190 km s" 1 


190-220 km i 


D n 4000 


0.0255 


0.0089 


0.0041 


0.0013 


0.0001 




0.0240 


0.0030 


-0.0020 


-0.0032 


-0.0020 


CNi 


0.0032 


0.0014 


0.0009 


0.0004 


0.0001 


CN 2 


0.0038 


0.0016 


0.0009 


0.0004 


0.0001 


Ca4227 


0.0058 


0.0059 


0.0050 


0.0035 


0.0015 


G4300 


0.0292 


0.0112 


0.0057 


0.0023 


0.0005 


H7F 


-0.0698 


-0.0255 


-0.0123 


-0.0045 


-0.0006 


Fe4383 


0.0765 


0.0456 


0.0324 


0.0195 


0.0074 


Ca4455 


0.0069 


0.0071 


0.0061 


0.0042 


0.0018 


Fc4531 


0.0298 


0.0153 


0.0100 


0.0055 


0.0019 


C 2 4668 


0.2019 


0.1160 


0.0807 


0.0478 


0.0178 


H/3 


0.0104 


0.0103 


0.0088 


0.0060 


0.0025 


Mgi 


0.0019 


0.0008 


0.0004 


0.0002 


0.0001 


Mg 2 


0.0054 


0.0023 


0.0013 


0.0006 


0.0002 


Mg 6 


0.0873 


0.0328 


0.0164 


0.0064 


0.0011 


Fc5270 


0.0474 


0.0303 


0.0222 


0.0137 


0.0053 


Fe5335 


0.0476 


0.0286 


0.0203 


0.0123 


0.0047 


Fc5406 


0.0382 


0.0187 


0.0118 


0.0064 


0.0021 


NaD 


0.1200 


0.0563 


0.0345 


0.0180 


0.0058 


TiOi 


-0.0017 


-0.0006 


-0.0002 


-0.0001 


0.0000 


TiQ 2 


0.0006 


0.0004 


0.0003 


0.0002 


0.0001 



Note. — Values given arc the scaled corrections for incompleteness, as described in §A. These arc subtracted 
from the measured index strengths given in Table 2. The same correction is applied to quiescent, weak, and 
strong LINER-likc emission composite spectra. All incompleteness corrections are given in the units of the 
corresponding index (CNi, CN2, Mgi, Mg2, TiOi, and Ti02 are in magnitudes, D.^4000 is dimcnsionless, and 
all others are in A). 



TABLE 7 
Abundance Trends with a 



Abundance = a + b log a 



Abundance 




b 








a 








[Fc/H] 







.451 


± 


0. 


034 


-1 


.14 


± 





.08 


[Mg/Fe] 







.362 


± 


0. 


.043 


-0 


.61 


± 


0. 


10 


[C/Fc] 







555 


± 


0. 


027 


-1 


06 


± 


0. 


06 


[N/Fc] 







.640 


± 


0. 


035 


-1 


,34 


± 


0. 


.08 


[Ca/Fe] 







130 


± 


0. 


.039 


-0. 


.29 


± 





09 


[Z/H] ([O/Fc] = 


= 0.0) 





.786 


± 


0. 


050 


-1 


.80 


± 





11 


[Z/H] ([O/Fe] = 


= 0.3) 





.705 


± 


0. 


051 


-1 


.47 


± 





11 


[Z/H] ([O/Fe] = 


= 0.5) 


0. 


.656 


± 


0. 


.052 


-1 


.23 


± 


0. 


12 


[Z/H] ([O/Fe] = 


= [Mg/Fe]) 


0. 


.916 


± 


0. 


.054 


-2. 


.00 


± 


0. 


12 


[Z/H] ([O/Fe] = 


= [Mg/Fe] 4 


- 0.3) 


.895 


± 


0. 


055 


-1 


.77 


± 


0. 


12 
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TABLE 8 

Trends with a: Comparison with Other Work 



log Age = a 


+ b log a 


This Work 


Trager et al. 


Thomas et al. 


Nelan et al. 


Smith et al. 


Bcrnardi et al. 


Abundance = 


a + b log a 




(2000b) 


(2005) 


(2005) 


(2007) 


(2003b) 


Age 


b 


0.35 ± 0.03 b 


0.6±0.2 C 


0.24 (0.32) d 


0.59 ±0.13 


0.52 ±0.06 


1.15 


[Fe/H] 


b 


0.45 ± 0.03 


0.48±0.12 c 








[Z/H] 


b 


0.79 ± 0.05° 


0.76±0.13 c 


0.55 (0.57) d 


0.53 ±0.08 


0.34 ±0.04 


0.38 


[a/Fe] a 


b 


0.36 ±0.04 


0.33 ±0.01 


0.28 (0.28) d 


0.31 ±0.06 


0.23 ±0.04 


0.32 


SSP Models' 




S06 


W94 


TMB03 


TMB03 


TMB03 


TMB03 



a Measured using Mg b in all works except Nelan et al.(2005) which uses Mg b and CNi as a sensitive indices. Corresponds to [Mg/Fc] 
in this work.* 5 From fit to H/3 age, quiescent galaxies only. c Trager et al. (2000b) posit the existence of an "age-metallicity hyperplanc" 
such that galaxies display a range of ages and [Z/H] (or [Fe/H]) at fixed a, with an anticorrelation between age and [Z/H]; they 
therefore give [Z/H] as a function of a and age. The age relation given here is from the Nelan et al.(2005) fit to the Trager et al. 
data, while the [Z/H] and [Fc/H] values include only the a portion of the effect on [Z/H] and [Fe/H] reported in Trager et al. (2005b), 
ignoring the effect of agc. d Thomas et al.(2005) report two sets of relations for galaxies in low and high density environments. Values 
for both environments arc show as: high density (low density). Calculated with [O/Fe] — 0.0, as described in §5.3. For other choices 
of [O/Fc], sec Table 5. f S06: Schiavon (2007), W94: Worthey et al. (1994), TMB03: Thomas, Maraston, & Bender (2003) 




FlG. 1. — (a) Ha and [On] EWs for all galaxies in the redshift range 0.06 < z < 0.08 and with both lines detected at > 2a. The 
distribution is bimodal (see Yan et al. 2006 for details). The EW(Ha)/EW([Oll]) criterion of equation 2 is used to define the LINER-like 
sample is shown as the solid line. Galaxies to the right of the solid line were discarded as "Low- [On] /Ha" contaminants with dusty star 
formation, Seyfert emission, or a mixture of Seyfert and star forming emission (Transition Objects). The dashed line shows the linear 
least squares fit of EW(Ha) as a function of EW([Oll]) which was used to convert a measured [On] EW into an Ha EW as part of the 
infill-correction process (see §2.2). (b) The same as (a) but with linear-scaled axes, (c) Color-magnitude diagram of all SDSS galaxies in 
the redshift range 0.06 < z < 0.08. The dashed line shows eq. 1 used to define the red sequence galaxies included in the sample, (d) 
Color-magnitude diagram of LINER-like and quiescent galaxies selected for the sample presented here. The dashed line shows the cut 
used to define the red sequence. The LINER-like and quiescent samples have similar color distributions (mostly red sequence galaxies) 
but different magnitude distributions. See text for details. In all panels, contour levels in all panels indicate the density of objects in the 
plot corresponding to, from light to dark, 3, 10, 30, 60, 100, 150, 200, 300, 400 galaxies per bin, where bin boundaries correspond to the 
tickmarks shown on the x and y axes. 
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Fig. 2. — BPT diagrams (Baldwin, Phillips, & Terlevich 1981) for galaxies with 0.06 < z < 0.08. Red points show those galaxies 
classified as "LINER-like" in our sample based on their EW([Oll])/EW(Ha) ratios; black points show all other galaxies. Only galaxies 
with all four emission lines detected at the 3cr level are plotted in each panel. The weak [Ol] emission line is undetected in many galaxies, 
thus (b) contains fewer galaxies than (a). Error bars in the upper left show the median error in line ratios separately for the LINER-like 
and other galaxies. The solid lines show the standard classification scheme: Seyferts have [Oni]/H/3 > 3 and [NllJ/Ha > 0.6, LINERs have 
[Oni]/H/3 < 3 and [Nil] /Ha > 0.6, and star forming galaxies have [Nll]/Has < 0.6. Dotted lines show the separation of Kewley et al. (2001) 
between star forming galaxies and active nuclei; the dashed line in (a) shows the Kauffmann et al. (2003) separation. The dash-dot line 
in (b) shows the Ho, Filippenko, & Sargent (1997) separation between LINERs and transition objects (TOs). Using the Kauffmann et 
al. (2003) separation between star forming galaxies and active nuclei, the standard separation between Seyferts and LINERs, and the Ho, 
Filippenko, & Sargent (1997) separation between LINERs and TOs, our LINER-like sample defined using EW([Oll])/EW(Ha) contains 
88.4% LINERs, 8.6% TOs, 1.9% Seyferts, and 1.1% star forming galaxies. 




EW([0II]) (A) 

Fig. 3. — Histogram of [On] EW in the quiescent galaxy sample. The unfilled histogram shows all galaxies that meet the selection 
criteria of §2.1. From these galaxies, a sub-sample of 2000 quiescent galaxies is chosen for comparison with the weak and strong LINER-like 
samples. The sub-sample, shown as the filled histogram, is constructed to have a symmetric gaussian distribution around EW([On]) = 0, 
and is truncated at ±3.12 A(the 2a spread of the total EW([Oll]) distribution) to limit contamination from undetected faint emission lines. 
See §2.3 for details. 
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Fig. 4. — Stacked spectra for the galaxy bins with 70 < a < 120 km s" 1 (the lowest a bins), with and without emission infill corrections. 
The black, green, and red lines show the quiescent, weak LINER-like, and strong LINER-like galaxy bins, respectively. Solid lines show 
the uncorrected spectra, while dotted lines show the emission infill correction. Differences are only visible in the Balmer lines where the 
corrections have been made. The correction to H<5 in particular is very small. Strong emission lines are indicated by vertical dashed lines 
and are labelled. Except for the emission lines, the spectra are remarkably similar. This illustrates the extremely high S/N of the stacked 
spectra and shows the wealth of information in the absorption features, as well as the subtlety of the spectral variations that are quantified 
in this paper. 



29 




-19 -20 -21 -SS -23 

-19 -20 -21 -22 -23-19 -20 -21 -22 -23-19 -20 -21 -22 -23 




-19 -20 -21 -22 -23-19 -20 -21 -22 -23-19 -20 -21 -22 -23 
Wfci T M* L i r 



Fig. 5. — Color-magnitude diagrams (CMDs) of quiescent and LINER-like red sequence galaxies. The top panel shows the CMD of all 
(quiescent and LINER-like) galaxies in the sample color-coded by aperture-corrected velocity dispersion (<r): <r = 70—120 km s" 1 in purple, 
a = 120-145 km s" 1 in blue, a = 145-165 km s _1 in green, a = 165-190 km s' 1 in gold, a = 190-220 km s _1 in salmon, and a = 220-300 
km s — 1 in red. The lower panels show similar CMDs with all galaxies in gray and overplotted colored contours for the distribution of the a 
range indicated. Colors correspond to the colors in the top panel. Below each contoured CMD is a histogram of the magnitude distribution 
of galaxies in that a bin. In all panels, the black dashed line shows the 100% completeness limit, while the solid black line shows the 
magnitude at which the sample is 50% incomplete. The histograms show that the highest a bins are effectively complete while the lower a 
bins are seriously incomplete. A correction for this effect is described in §A. Higher a galaxies are more luminous and redder than lower a 
galaxies, but there is a large spread in both magnitude and color at fixed <r, typically larger than the separation between a bins. At fixed 
a, there is no clear color-magnitude relation. 
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Fig. 6. — Lick index measurements in the stacked spectra as a function of the mean a and emission line properties of the galaxies in each 
stacked spectrum. Quiescent galaxies are shown in black, weak LINER-like galaxies in green, and strong LINER-like galaxies in red. Error 
bars show the \a line strength measurement errors. Solid lines show linear least-squares fits of index strength onto a, independently for 
each bin in emission line properties, (a) Indices strongly sensitive to age. (b) Indices dominated by iron absorption lines. In each panel, 
the left column shows raw measured indices, while the right column shows the indices corrected for incompleteness in the lower a bins (see 
§A). The Balmer line strengths decrease with increasing a while all other line strengths increase with increasing cr. This is consistent with 
higher <r galaxies being older and/or more metal rich than lower <r galaxies. Galaxies with LINER-like emission have stronger Balmer lines 
and weaker metal lines, consistent with LINER-like galaxies having younger age and/or lower metallicities than quiescent galaxies, (c) 
and (d) Indices strongly affected by non-solar abundance ratios. The same trends are visible as in panel (b), that higher a galaxies have 
stronger absorption lines than lower a galaxies while LINER-like galaxies have weaker absorption lines than quiescent galaxies (with the 
exception of Na D, where LINER-like galaxies have stronger absorption than quiescent galaxies; see §3.2 for details). 




cont. 
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Fig. 7. — The effect of emission infill and various corrections for emission infill on Balmer absorption line strengths. As before, quiescent 
galaxies are shown in black, weak LINER-like galaxies in green, and strong LINER-like galaxies in red; error bars show measurement errors 
for the stacked spectra, and solid lines are linear least squares fits of line strength onto <r. (a-c) Balmer absorption line measurements with 
no infill correction applied. Note that in (c) H/3 is substantially filled in by emission in the LINER-like galaxies while in (a) H<5p is only 
marginally affected. Panels (d-m) show the effects of making a range of reasonable emission infill corrections, (d-f) show the Balmer lines 
when the Balmer decrement measured by Yan et al. (2006) is used. H/3 line strengths vary little between quiescent and LINER-like galaxies, 
while H5p line strengths vary substantially, (g— i) show the Balmer lines using the Balmer decrement measured in red sequence LINERs of 
Ho, Filippenko, & Sargent (1997). These two choices of the Balmer decrement likely bracket the true mean value of the Balmer decrement 
in these systems, (j— m) show the Balmer lines using a "maximal" Balmer decrement, with Ha/H/3 characteristic of stellar-ionized Hn 
regions and no reddening. The trend of stronger Balmer absorption in LINER-like galaxies exists for all reasonable values of the Balmer 
decrement, and exists even in the uncorrected H8p, where emission infill is minimal. 
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Fig. 8. — Comparison of Lick index measurements for 3 stacked spectra with a simple stellar population model from Schiavon (2007). 
The triangle, star, and square, show measurements for the 70 < a < 120 km s — 1 , 145 < a < 165 km s — 1 , and 220 < a < 300 km s~ 1 
quiescent galaxy bins, respectively. Solid lines connect models of the same [Fe/H], from left to right: —1.3, —0.7, —0.4, 0.0, and +0.2. 
Dotted lines connect models of the same age, with ages from top to bottom: 1.2, 2.2, 3.5, 7.9, and 14.1 Gyr, as labelled in the top left 
plot. For reference, the [Fe/H] = 0.0 and age = 7.9 Gyr lines are shown as the long-dash and short-dash lines, respectively. The model 
has abundance ratios [Mg/Fe] = [Na/Fe] = [Si/Fc] = [Ti/Fe] = 0.120, [C/Fc] = 0.100, [N/Fe] = -0.002, and [Ca/Fe] = 0.010, calculated 
with a solar isochrone and with [O/Fc] = [Cr/Fc] = 0.0. This is the best-fitting model for the 70 < a < 120 km s _1 data (triangle), as 
determined by EZ_Ages (Graves & Schiavon 2007). Note that the data points for the 70 < a < 120 km s _1 stacked spectrum fall at the 
same values of age and [Fe/H] in each panel, as expected since this model gives consistent parameters for all indices. The data points from 
the 145 < a < 165 km s _1 (star) and 220 < a < 300 km s _1 (square) stacked spectra do not fall at the same values in each panel; higher 
values of [Mg/Fc], [C/Fe], and [N/Fe] are needed to make the age and [Fe/H] measured in Mg 6, C24668, and CN2 match the other panels, 
suggesting that the higher a galaxies are enhanced in [Mg/Fe], [C/Fe], and [N/Fe] over the 70 < a < 120 km s _1 galaxies. 
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Fig. 9. — Index-index plots for indices not used by EZ_Ages (Graves & Schiavon 2007) in the abundance fitting process. All lines and 
symbols are as in Figure 8. The age and [Fc/H] values for the 70 < a < 120 km s _1 galaxies (triangle) are in good agreement with those 
seen in Figure 8, with the exception of G4300 and Yi-fp. Also, the higher <r galaxies have stronger Mg2 and CNi than predicted by this 
model, suggesting they have higher [Mg/Fe], [C/Fe], and [N/Fe] than the 70 < a < 120 km s _1 galaxies, as was also implied by Figure 8. 
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Fig. 10. — Same as Figure 8 except that grids are computed using the best-fit abundances for the quiescent 220 < a < 300 km s 1 
data, as determined by EZ_Ages. This model has [Mg/Fc] = [Na/Fc] = [Si/Fe] = [Ti/Fe] = 0.260, [C/Fe] = 0.270, [N/Fc] = 0.220, and 
[Ca/Fe] = 0.040, calculated with a solar isochrone and with [O/Fe] = [Cr/Fe] = 0.0. The triangle, star, and square show the data for the 
70 < a < 120 km s —1 , 145 < a < 165 km s — 1 , and 220 < cr < 300 km s _1 quiescent stacked spectra, respectively. Here, note that the data 
point for the 220 < <r < 300 km s _1 stacked spectrum falls in the same position on the grid in each index-index plot, indicating that this 
is the correct abundance pattern for this data. The lower a galaxies, by contrast, are clearly less enriched in Mg, CH, and CN than this 
model as shown by their location in the plots of H/3 against Mg b, C24668, and CN2. 
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Fig. 11. — Same as Figure 9 except that grids are computed using the best-fit abundances for the quiescent 220 < a < 300 km s~ 1 
data, showing indices not used in the fitting process. The plots of H/3 against Fe4383, Mg2, and CNi give ages and [Fe/H] values for the 
220 < a < 300 km s" 1 data which are consistent with those in Figure 10, indicating that this abundance pattern reproduces even indices 
which are not used in the fitting process. The G4300 index is inconsistent with other indices, as discussed in the text. The ages measured 
from and H<5p are slightly younger than the ages measured from H/3. 
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Fig. 12. — SSP ages as a function of a and emission line properties. As in previous plots, quiescent, weak, and strong LINER-likc 
galaxies are shown in black, green, and red, respectively. Error bars show the uncertainty in the ages due to measurement errors in the 
Balmer lines and (Fe) index only. Solid lines show linear least squares fits of age onto a for each data set. Ages are determined separately 
for the three Balmer lines H/3, H7J?, and H<5p. The left and center columns show the ages measured with the Ho, Filippenko, & Sargent 
(1997, our preferred correction) and Yan et al. (2006) Balmer decrements, respectively, to illustrate the effect of different infill corrections. 
Both are determined by EZ_Ages (Graves & Schiavon 2007) using a solar-scale isochrone. The right column shows the effect on the age 
measurements when an a-enhanced isochrone and [O/Fe] = +0.3 are adopted. Higher <r galaxies have older ages than lower cr galaxies, 
and LINER-like galaxies are younger than quiescent galaxies at fixed <r, independent of the choice of Balmer decrement, stellar population 
model isochrone, or Balmer line used to make the measurement. However age differences are more pronounced using higher-order Balmer 
lines, as expected if the stellar populations have a distribution of ages. A smaller Balmer decrement (i.e., larger infill correction) leads to 
larger differences between LINER-like and quiescent ages. The use of an a-enhanced isochrone produces systematically younger ages in 
H/3, but qualitatively similar trends with <r and emission line properties. 



8<S 



Graves et al. 



1.10 r 

1.05 r 




- LOO: 

U 0.95 : 

S 0.90 - 

0.85 : 
0.80 L 



Fig. 13. — Median 0,1 (g — r) colors as a function of mean cr for the galaxies that went into each stacked spectrum. The quiescent 
galaxies are shown in black, the weak LINER-likc galaxies in green, and the strong LINER-like galaxies in red. Solid lines show a linear 
least-squares fit of color onto log a. The error bars show the la spread in the color values for each bin of galaxies. Although the spectral 
data show that the LINER-like galaxies are typically younger than their quiescent counterparts, the LINER-like galaxies on average have 
slightly redder colors. This is not a metallicity effect (see Figure 14) but is probably an indication that the LINER-like galaxies are typically 
more dust-reddened than the quiescent galaxies, consistent with the larger amounts of interstellar NaD absorption seen in the LINER-like 
galaxies in Figure 6. See §5.1 for details. 
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Fig. 14. — (a-e) [Fe/H] and abundance ratios for quiescent and LINER-like galaxies as a function of a. Quiescent, weak, and strong LINER- 
like galaxies are shown in black, green, and red, respectively. Error bars show the uncertainties in each abundance due to measurement 
errors in the Lick indices as described in §5.3. The dotted line in each panel shows solar abundance. There are no significant variations in 
abundances or abundance ratios between quiescent and LINER-like galaxies. The solid line show a linear least squares fit of abundance onto 
a for quiescent and LINER-like galaxies simulatenously. Abundances are computed using a solar-scale isochrone and assuming [O/Fe] = 
0.0, with Na, Si, and Ti set to track Mg, and [Cr/Fe] = 0.0. The abundances computed using the o-enhanced isochrone are not significantly 
different. The dashed line shows the effect of computing abundances with [O/Fe] = +0.3. All abundances increase with cr, [Fe/H], [C/Fc], 
and [N/Fc] most strongly. All galaxies are Mg-enhanced. (f) Linear least squares fits of [Z/H] as a function of a. The black solid, dashed, 
and dash-dot lines show [Z/H] computed assuming [O/Fe] = 0.0, [O/Fe] = +0.3, and [O/Fe] = +0.5, respectively. The blue solid and 
dashed lines show [Z/H] computed assuming [O/Fe] = [Mg/Fe], and [O/Fe] = [Mg/Fe] + 0.3, respectively. Total metallicity is dominated 
by [O/Fe] which cannot be measured and is thus highly uncertain. 
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Fig. 15. — Lick index line strengths in stacked spectra of quiescent red sequence galaxies, as a function of a and magnitude. Each 
bin in a has been divided into faint, mid-magnitude, and bright sub-bins, with equal numbers of galaxies in each sub-bin. The faintest 
galaxies at each a are shown in blue, the mid-magnitude galaxies are shown in green, and the brightest galaxies are shown in red. Error 
bars show the la index measurement errors for each stacked spectrum. Solid lines show linear least-squares fits of index strength onto a, 
independently for each bin in magnitude. As in Figure 6, panel (a) shows strongly age-sensitive indices, panel (b) shows indices dominated 
by iron absorption lines, and panels (c) and (d) show indices that are sensitive to non-solar abundance ratios. Faint galaxies appear to 
have similar ages but lower metal abundances when compared to mid-magnitude and bright galaxies at the same a. 
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Fig. 16. — [Fe/H], abundance ratios, and ages of stacked spectra, as a function of galaxy absolute magnitude. Black stars are the 
luminosity-binned data from appendix B. Only the stacked spectra for quiescent galaxies arc shown. The solid black line shows a linear 
least-squares fit of the stellar population parameter onto magnitude for the luminosity-binned data. Grey triangles show the results from 
Schiavon (2007) for the Eiscnstein et al. (2003) stacked spectra. The dotted grey line shows a linear least-squares fit for the Eisenstein et al. 
(2003) stacked spectra. Solar abundance and abundance patterns arc shown for reference as the dashed line in each panel. The Eiscnstein 
ct al. (2003) data (triangles) show younger ages, slightly higher [Fc/H], and slightly steeper abundance relations versus magnitude than do 
the luminosity-binned spectra from our sample (stars). 
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Fig. 17. — Selected Lick index line strengths as a function of galaxy magnitude. Black stars are the luminosity-binned data from appendix 
B. Only the stacked spectra for quiescent galaxies are shown. Grey triangles show indices measured in the Eisenstcin ct al. (2003) stacked 
spectra. In both samples, Balmer lines weaken as galaxy luminosity increases with the exception of H/3, which behaves similarly to }18p 
and in our data but is constant in the Eisenstcin ct al. (2003) data. H<5^ and }l-fp are stronger in the Eisenstein et al. (2003) data 

than in ours, consistent with younger mean ages. In contrast to the Balmer lines, the other Lick indices are very similar between the two 
data sets. The metal lines may be slightly weaker in the Eisenstein et al. (2003) data than in ours, again consistent with younger ages. 
The consistency of the metal absorption line strengths between the two samples implies that the differences in the measured abundance 
trends seen in Figure 16 are caused by the difference in measured Balmer line strengths, rather than intrinsic differences in mctallicity. 



